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Abstract 

Galaxy redshift surveys have achieved significant progress over the 
last couple of decades. Those surveys tell us in the most straightforward 
way what our local universe looks like. While the galaxy distribution 
traces the bright side of the universe, detailed quantitative analyses of 
the data have even revealed the dark side of the universe dominated by 
non-baryonic dark matter as well as more mysterious dark energy (or 
Einstein's cosmological constant). We describe several methodologies of 
using galaxy redshift surveys as cosmological probes, and then summarize 
the recent results from the existing surveys. Finally we present our views 
on the future of redshift surveys in the era of Precision Cosmology. 
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1 Introduction 



Nowadays the exploration of the universe can be performed by a variety of obser- 
vational probes and methods over a wide range of the wavelengths; temperature 
anisotropy map of the cosmic microwave background (CMB), the Hubble dia- 
grams (i.e., redshift-magnitude relations) of nearby galaxies and distant Type 
la supernovae, wide-field photometric and spectroscopic surveys of galaxies, 
power spectrum and abundances of galaxy clusters in optical and X-ray bands 
combined with the radio observation through the Sunyaev-Zel'dovich effect, 
deep surveys of galaxies in submm, infrared, and optical bands, quasar surveys 
in radio and optical, strong and weak lensing of distant galaxies and quasars, 
high-energy cosmic rays and so on. Undoubtedly gamma-rays, neutrinos, grav- 
itational radiations will join the above already crowded list. 

Among those, optical galaxy redshift surveys are the most classical. Indeed 
one may phrase that the modern observational cosmology started with a sort of 
a galaxy redshift survey by Edwin Hubble. Still galaxy redshift surveys are of 
vital importance in cosmology in the 21st century for various reasons: 

Redshift surveys have unprecedented quantity and quality: 

The mimbcrs of galaxies and quasars in the spectroscopic sample of Two 
Degree Field (2dF: see §6.1.1) are ~ 250, 000 and - 30, 000, and will reach 
~ 800, 000 and 100, 000 upon completion of the on-going Sloan Digital Sky 
Survey (SDSS: see §6.1.2). These unprecedented munbers of the objects 
as well as the homogeneous selection criteria enable the precise statistical 
analysis of their distribution. 

The universe at ^; ~ 1000 is well specified: 

The first-year WMAP ( Wilkinson Microwave Anisotropy Probe) data [6] 
among others have established a set of cosmological parameters. This 
may be taken as the initial condition of the universe from the point-of- 
vicw of the structure evolution toward z = 0. In a sense, the origin of 
the universe at z ~ 1000 and the evolution of the universe after the epoch 
are now equally important, but well separable questions that particle and 
observational cosmologists focus on, respectively. 

Gravitational growth of dark matter component is well understood: 

In addition, extensive numerical simulations of structure formation in the 
universe has significantly advanced our understanding of the gravitational 
evolution of dark matter component in the standard gravitational insta- 
bility picture. In fact, we even have very accurate and useful analytic 
formulae to describe the evolution deep in its nonlinear regime. Thus we 
can now directly address the evolution of visible objects from the analy- 
sis of their redshift surveys separately from the nonlinear growth of the 
underlying dark matter gravitational potentials. 
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Formation and evolution of galaxies: 

In the era of precision cosmology among others, scientific goals of re- 
searches out of galaxy redshift surveys are gradually shifting from inferring 
a set of values of cosmological parameters using galaxy as their probes to 
understanding the origin and evolution of galaxy distribution given a set 
of parameters accurately determined by the other probes like CMB and 
supernovae. 

With the above in mind, we will attempt to summarize what wc have learned 
so far from galaxy redshift surveys and then to describe what will be done 
with future data. The review is organized as follows. We first present a brief 
overview of the Fricdmann model and gravitational instability theory in §2. 
Then we describe the non-Gaussian nature of density fluctuations generated by 
the nonlinear gravitational evolution of the primordial Gaussian field (§3). Next 
we disciiss the spatial biasing of galaxies relative to the underlying dark matter 
distribution in §4. Our understanding of biasing is still far from complete, and 
its description is necessarily empirical and very approximate. Nevertheless this 
is one of the most important ingredients for proper interpretation of galaxy 
redshift surveys. Section 5 introduces general relativistic effects which become 
important especially for galaxies at high redshifts. We present the latest results 
from the two currently largest galaxy redshift surveys, 2dF (Two Degree Field) 
and SDSS (Sloan Digital Sky Survey), in §6. Finally §7 is devoted to summary 
of the present knowledge of our universe and our personal view of the future 
direction of cosmological researches in the new millennium. 
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2 Clustering in the expanding universe 



2.1 The cosmological principle 

Our current universe exhibits a wealth of nonhnear structures, but the zero-th 
order description of our universe is based on the assumption that the universe is 
homogeneous and isotropic smoothed over sufficiently large scales. This state- 
ment is usually referred to as the cosm,ological principle. In fact, the cosmologi- 
cal principle was first adopted when observational cosmology was in its infancy; 
it was then little more than a conjecture, embodying 'Occam's razor' for the 
simplest possible model. 

Rudnicki (1995) [73] summarized various forms of cosmological principles in 
modern-day language, which were stated over different periods in human history 
based on philosophical and aesthetic considerations rather than on fundamental 
physical laws. 

(i) The ancient Indian cosmological principle: 

The Universe is infinite in space and tirn,e and is infinitely heterogeneous. 

(ii) The ancient Greek cosmological principle: 

Our Earth is the natural center of the Universe. 

(ill) The Copernican cosmological principle: 

The Universe as observed from any planet looks much the same. 

(iv) The (generalized) cosmological principle: 

The Universe is (roughly) hom.ogeneous and isotropic. 

(v) The perfect cosmological principle: 

The Universe is (roughly) homogeneous in space and time, and is isotropic 
in space. 

(vi) The anthropic principle: 

A human being, as he/she is, can exist only in the Universe as it is. 

We note that the ancient Indian principle may be viewed as a sort of 'fractal 
model'. The perfect cosmological principle led to the steady state model, which 
although more symmetric than the (generalized) cosmological principle, was 
rejected on observational grounds. The anthropic principle is becoming popular 
again, e.g. in 'explaining' the non-zero value of the cosmological constant. 

Like with any other idea about the physical world, we cannot prove a model, 
but only falsify it. Proving the homogeneity of the Universe is in particular dif- 
ficult as we observe the Universe from one point in space, and we can only 
deduce isotropy indirectly. The practical methodology we adopt is to assume 
homogeneity and to assess the level of fiuctuations relative to the mean, and 
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hence to test for consistency with the underlying hypothesis. If the assumption 
of homogeneity turns out to be wrong, then there are numerous possibiUties for 
inhomogeneous models, and each of them must be tested against the observa- 
tions. 

For that purpose, one needs observational data with good quality and quan- 
tity extending up to high redshifts. Let us mention some of those. 

The CMB Fluctuations Ehlers, Garcn and Sachs (1968) [20] showed that 
by combining the CMB isotropy with the Copernican principle one can 
deduce homogeneity. More formally their theorem (based on the Liouville 
theorem) states that "If the fundamental observers in a dust space-time 
see an isotropic radiation field, then the space-time is locally given by 
the Friedman-Robertson- Walker (FRW) metric". The COBE (COsmic 
Background Explorer) measurements of temperature fluctuations AT/T = 
10-5 on scales of 10° give via the Sachs Wolfe effect (AT/T = ^A^/c^) 
and Poisson equation rms density fluctuations of ^ ^ 10"^ on 1000 Mpc 
(e.g. [96]), which implies that the deviations from a smooth Universe are 
tiny. 

Galaxy Redshift Surveys The distribution of galaxies in local rcdshift sur- 
veys is highly clumpy, with the Supergalactic Plane seen in full glory. 
However, deeper surveys like 2dF and SDSS (see §6) show that the fluctu- 
ations decline as the length-scales increase. Peebles (1993) [69] has shown 
that the angular correlation functions for the Lick and APM (Automatic 
Plate Measuring) surveys scale with magnitude as expected in a universe 
which approaches homogeneity on large scales. While rcdshift surveys can 
provide interesting estimates of the fluctuations on intermediate scales 
(e.g. [72]), the problems of biasing, evolution and iv'-correction, would 
limit the ability of those redshift surveys to 'prove' the Cosmological Prin- 
ciple. Despite these worries the measurement of the power spectrum of 
galaxies derived on the assumption of an underlying FRW metric shows 
good agreement with the A-CDM (cold dark matter) model. 

Radio Sources Radio sources in siirveys have typical median redshift z ^ 1, 
and hence are useful probes of clustering at high redshift. Unfortunately, 
it is diflicult to obtain distance information from these surveys: the radio 
luminosity function is very broad, and it is difficult to measure optical 
redshifts of distant radio sources. Earlier studies claimed that the distri- 
bution of radio sources supports the 'Cosmological Principle'. However, 
the wide range in intrinsic luminosities of radio sources would dilute any 
clustering when projected on the sky. Recent analyses of new deep ra- 
dio surveys suggest that radio sources are actually clustered at least as 
strongly as local optical galaxies. Nevertheless, on the very large scales 
the distribution of radio sources seems nearly isotropic. 
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The XRB The X-ray Background (XRB) is hkely to be due to sources at high 
redshift. The XRB sources are probably located at redshift z < 5, making 
them convenient tracers of the mass distribution on scales intermediate 
between those in the CMB as probed by COBE, and those probed by op- 
tical and IRAS redshift surveys. The interpretation of the results depends 
somewhat on the nature of the X-ray sources and their evohition. By 
comparing the predicted multipoles to those observed by HEAOl, Scharf 
et al. (2000) [75] estimate the amplitude of fluctuations for an assumed 
shape of the density fluctuations. The observed fluctuations in the XRB 
are roughly as expected from interpolating between the local galaxy sur- 
veys and the COBE and other CMB experiments. The rms fluctuations 
^ on a scale of ~ eOO/i-^Mpc are less than 0.2 %. 

Since the (generalized) cosmological principle is now well supported by the 
above observations, we shall assume below that it holds over scales / > 100/i~^ 
Mpc. 

The rest of the current section is devoted to a brief review of the homoge- 
neous and isotropic cosmological model. Further details may be easily found in 
standard cosmology textbooks [93, 62, 69, 64, 10, 63]. 

The cosmological principle is mathematically paraphrased as that the metric 
of the universe (in its zero-th order approximation) is given by 

ds^ = -de + a{tf ( + ^''^^^ + xHmHd<j)^^ , (1) 

where x is the comoving coordinate, and we use the units in which the light 
velocity c = 1. The above Robertson - Walker metric is specified by a constant 
K, the spatial curvature, and a function of time a{t), the scale factor. 

The homogeneous and isotropic assumption also implies that the energy 
momentum tensor of the matter field, T^^, should take the form of the ideal 
fiuid: 

T/^i^ = iP + P)u^u^ +pg^u, (2) 

where is the 4-velocity of the matter, p is the mean energy density, and p is 
the mean pressure. 

2.2 From the Einstein equation to the Friedmann equation 

The next task is to write down the Einstein equation: 

Rfii^ - ]^Rg^u + ^Qui^ = SnGTf^,, (3) 

using equations (1) and (2). In this case one is left with the following two 
independent equations: 



1 (Pa A-kG , „ , A 



for the three independent functions, a{t), p{t), and p{t). 
Differentiation of equation (4) with respect to t yields 

a AnG / . a „ \ A 



a 3 V d / 3 
Then ehminating a with equation (5), one obtains 

P=-^l{p+p). (7) 

This can be easily interpreted as the 1st law of the thermodynamics dQ = 
dU — pdV = d{pa?) — pd{a^) = in the present context. Equations (4) and 
(7) arc often used as the two independent basic equations for a{t), instead of 
equations (4) and (5). 

In either case, however, one needs another independent equation to solve for 
a{t). This is usually given by the equation of state of the form p = p{p). In 
cosmology, the following simple relation: 

p = wp. (8) 

While the value of w may change as redshift in principle, it is often assumed that 
w is independent of time just for simplicity. Then substituting this equation of 
state into equation (7) immediately yields 

pcxa-3(i+"'\ (9) 

The non-relativistic matter (or dust), ultra-relativistic matter (or radiation), 
and the cosmological constant correspond to w = 0, 1/3, and —1, respectively. 

If the universe consists of different fluid species with Wi (i = 1,...,N), 
equation (9) still holds independently as long as they do not interact with each 
other. If one denotes the present energy density of the i-th component by pi^o, 
then the total energy density of the universe at the epoch corresponding to the 
scale factor of a{t) is given by 

where the present value of the scale factor, ag, is set to be unity without loss of 
generality. Thus, equation (4) becomes 

IdaX^ _ 8'kG -^ Pi,o K A 

adt) ~ 3 ^ a3(i+«'0 a2 + 3 ■ ^ ' 
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Note that those components with Wi = —1 may be equivalent to the conventional 
cosmological constant A at this level, although they may exhibit spatial variation 
unlike A. 

Evaluating equation (11) at the present epoch, one finds 

SttG v-^ A 
^0 =— E^».°-^+ 3' (12) 

i=l 

where Hq is the Hubble constant at the present epoch. The above equation is 
usually rewritten as follows: 



K = H^Qk = Hi + 0^-1^ 



(13) 



where the density parameter for the i-th component is defined as 



^'.0 = -^JJ2 Pi'"' (14) 



and similarly the dimensionless cosmological constant is 

A 



^^A ^ (15) 



3Hi 



Incidentally equation (13) clearly illustrates the Mach principle in the sense that 
the space curvature is simply determined by the amount of matter components 
in the universe. In particular, the flat universe {K = 0) implies that the sum of 
the density parameters is unity: 



N 

Y^ni,, + n^ = i. (16) 

i=l 

Finally the cosmic expansion is described by 

I) ' = ^0 f E ^ + 1 - E - + ^^-') ■ (17) 

^ \«— 1 i—l / 

As will be shown below, the present universe is supposed to be dominated by 
non-relativistic matter (baryons and collisionless dark matter) and the cosmo- 
logical constant. So in the present review, we approximate equation (17) as 



dtj H^l^ + i-n^-n^ + n^a'] (is) 



unless otherwise stated. 
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2.3 Expansion law and age of the universe 

Equation (18) has the following analytic solutions in several simple but practi- 
cally important cases. 



(a) Einstein - de Sitter model (0„ = 1, Qa = 0) 

2/3 



a{t) = { A 



in 



2 



(19) 



(b) Open model with vanishing cosmological constant (r2„ < l,f2A = 0) 

fi^ „ a„ 



Hoto 



2(1 - n„) 
1 



(cosh^-l), Hot 
m 



2(1 - ri„)3/2 

2 - n„ + 2V1 - 0„ 



1 - 2(1 - f2„)3/2 



(sinh6l- 6») (20) 
(21) 



(c) Spatially-flat model with cosmological constant < 1, = 1 — f^m) 



a{t) = 

Hoto = 



1 - n„ 
1 



1/3 r 



sinh —Hot 



2 - 0„ + 2V1 - Ct„ 



2/3 



3V1 - f^m 

In the above, to denotes the present age of the universe: 

X dx 



t = — c 

Ho Jo 



^n^x + (1 - f2„ - nA)a:2 + f2Aa;4 



(22) 
(23) 

(24) 



and we adopt the initial condition that a = at < = 0. The expression clearly 
indicates that to increases as Q„ decreases and/or f^A increases. Figure 1 plots 
the scale factor as a function of Ho{t — to), and Table 1 summarizes the age of 
the universe. 



Table 1: The present age of the universe [(/i/0.7) ^ Gyr] 





Open model [Q.^ = 0) 


Spatially-flat model (Qa = 1 — f^m) 


1.0 


9.3 


9.3 


0.5 


10.5 


11.6 


0.3 


11.3 


13.5 


0.1 


12.5 


17.8 


0.01 


13.9 


28.0 
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Figure 1: Evolution of the cosmic scale factor as a function of Ho{t — t,,). The 
present value of the scale factor Oq is set to unity; solid line: (0„,riA) = 
(0.3,1.7), dotted Hne: (0.3,0.7), dot-dashed line: (0.3,0.0), long dashed Hne: 
(1.0,0.0), short dashed line: (3.0,0.0). 



2.4 Einstein's static model and Lemaitre's model 

So far we have shown that solutions of the Einstein equation are dynamical in 
general, i.e., the scale factor a is time-dependent. As a digression, let us ex- 
amine why Einstein once introduced the A-term to obtain a static cosmological 
solution. This is mainly important for historical reasons, but also is interesting 
to observe how the operationally identical parameter (the A-term, the cosmo- 
logical constant, the vaccum energy, the dark energy) shows up in completely 
different contexts in the course of the development of cosmological physics. 

Consider first the case of A = in equations (4) and (5). Clearly the 
necessary and sufficient condition that the equations admit the solution of a = 
const, is given by 

P^-^P-^^- (25) 

Namely, any static model requires that the universe is dominated by matter with 
either negative pressure or negative density. This is physically unacceptable as 
long as one considers normal matter in the standard model of particle physics. 
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If A ^ on the other hand, the condition for the static solution is 



K = 4TTG{p + p)a^, A = AwG{p + 3p), 



(26) 



yielding 



1 



3K 



A 



SttG V a 

Thus both p and p can be positive if 

K 



SttG 



A 



K 



<A< 



3^ 



In particular, if p = 0, 



A 

4^' 



(27) 



(28) 



(29) 



This represents the closed universe (with positive spatial curvature), and corre- 
sponds to Einstein's static model. 

The above static model is a special case of Lemaitre's universe model with 
A > and K > 0. For simplicity, let us assume that the universe is dominated 
by non-relativistic matter with negligible pressure, and consider the behavior 
of Lemaitre's model. First define the values of the density and the scale factor 
corresponding to Einstein's static model: 



47rG' 



K = Aal. 



(30) 



In order to study the stability of the model around the static model, consider a 
model in which the density at a = is a factor of a(> 1) larger than pe. Then 



apE I — 

a 



aA 
4^ 



and equations (4) and (5) reduce to 

.2 A / 2aa% 



a 
a 



3" \a 
A 
3" 



3a J + a' 

3' 



(31) 

(32) 
(33) 



For the period of a ^ Ce, equation (32) indicates that a cx t^/^ and the 
universe is decelerating (a < 0). When a reaches a^/^ag, takes the mini- 
mum value Aa|(a^/^ — 1) and the universe becomes accelerating (a > 0). Fi- 
nally the universe approaches the exponential expansion or de Sitter model: 
a oc exp{t-^/ A/3). If a becomes closer to unity, the minimum value reaches zero 
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and the expansion of the universe is effectively frozen. This phase is called the 
coasting period, and the case with a = 1 corresponds to Einstein's static model 
in which the coasting period continues forever. A similar consideration for a < 1 
indicates that the universe starts collapsing (d^ = 0) before a = 0. Thus the 
behavior of Lemaitre's model is crucially different if a is larger or smaller than 
unity. This suggests that Einstein's static model (a = 1) is unstable. 



2.5 Vacuum energy as an effective cosmological constant 

So far we discussed the cosmological constant introduced in the l.h.s. of the 
Einstein equation. Formally one can move the A-term to the r.h.s. by assigning 

This effective matter field, however, should satisfies an equation of state of 
p = —p. Actually the following example presents a specific example for an 
effective cosmological constant. Consider a real scalar field whose Lagrangian 
density is given by 

-C=^5,..a'^W-n<^)- (35) 
Its energy-momentum tensor is 

Ti^iy = dti4>dv4> - i^9nv, (36) 
and if the field is spatially homogeneous, its energy density and pressure are 

P0 = \^'' + y{4'l P4> = \^'' - V{4>). (37) 

Clearly if the evolution of the field is negligible, i.e., (j)^ <^ w —p^ and 

the field acts as a cosmological constant. Of course this model is one of the 
simplest examples, and one may play with much more complicated models if 
needed. 

If the A-term is introduced in the l.h.s., its should be constant to satisfy 
the energy-momentum conservation T^^^ = 0. Once it is regarded as a sort of 
matter field in the r.h.s, however, it does not have to be constant. In fact, the 
above example shows that the equation of state for the field has w = —1 only 
in special cases. This is why recent literature refers to the field as dark energy 
instead of the cosmological constant. 



2.6 Gravitational instability 

We have presented the zero-th order description of the universe neglecting the 
inhomogeneity or spatial variation of matter inside. Now we are in a position 
to consider the evolution of matter in the universe. For simplicity we focus on 
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the non-relativistic regime where Newtonian approximation is vahd. Then the 
basic equations for the self-gravitating fluid are given by the continuity equation, 
Euler's equation and the Poisson equation: 

^+V-{pu) = 0, (38) 

du 1 

— + {u-V)u = --Vj3-V$, (39) 

= AttGp. (40) 

Wo would like to rewrite those equations in the comoving frame. For this pur- 
pose, we introduce the position x in the comoving coordinate, the peculiar 
velocity, v, density fluctuations S{t,x), and the gravitational potential (l>{t,x) 
which are deflned respectively as 

« = if)' 

V = a{t)x, (42) 
- ^-1. (43) 
(f){t,x) = ^ + ^aa\xf. (44) 
Then equations (38) to (40) reduce to 

5 + -V ■ [{1 + 6)v] = 0, (45) 
a 

V + -(v-'V)v + -v = — —Vp--V(j), (46) 
a a pa a 

-- AirGpa^S, (47) 



72, 



where the dot and V in the above equations are the time derivative for a given 
X and the spatial derivative with respect to x, i.e., deflned in the comoving 
coordinate (while those in cqs.[38] to [40] arc defined in the proper coordinate). 

A standard picture of the cosmic structure formation assumes that the ini- 
tially tiny amplitude of density fluctuation grow according to equations (45) to 
(47). Also the universe smoothed over large scales approaches a homogeneous 
model. Thus at early epochs and/or on large scales, the nonlinear effect is small 
and one can linearize those equations with respect to 6 and v: 

S + iv-v = 0, (48) 
a 

V + -v = -^VS---V(j>, (49) 
a a a 

= ^nGpa'^6, (50) 
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where = [dp /dp) is the sound velocity. 

As usual, we transform the above equations in k space using 

5j^{t) = — j 6{t,x) exp{ik ■ x) dx. (51) 

Then the equation for Sj^ reduces to 



If the signature of the third term is positive, 5^ has an unstable, or, monotoni- 
cally increasing solution. This condition is equivalent to the Jeans criterion: 

(53) 

namely, the wavelength of the fluctuation is larger than the Jeans length Aj 
which characterizes the scale that the sound wave can propagate within the dy- 
namical time of the fluctuation yJ-K/Gp. Below the scale the pressure wave can 
suppress the gravitational instability, and the fluctuation amplitude oscillates. 




2.7 Linear growth rate of the density fluctuation 

Most likely our universe is dominated by coUisionless dark matter, and thus Aj 
is negligibly small. Thus at most scales of cosmological interest, equation (53) 
is well approximated as 

h + - 47rGp5^ = 0. (54) 

For a given set of cosmological parameters, one can solve the above equation by 
substituting the expansion law for a{t) as described in §2.3. Since equation (54) 
is the second-order differential equation with respect to t, there are two inde- 
pendent solutions; a decaying mode and a growing mode which monotonically 
decreases and increases as t, respectively. The former mode; bc;c,omes negligibly 
small as the universe expands, and thus one is usually interested in the growing 
mode alone. 

More specifically those solutions are explicitly obtained as follows. First note 
that the l.h.s. of equation (18) is the Hubble parameter at t, H{t) = d/a: 

= ^f^„(l + z)3 + (l-0„-f]J(l + z)2 + i7^] . (56) 
The flrst and second differentiation of equation (55) with respect to t yields 

,/ n 1-n 

2HH = Hi (-3^ - 2 H. (57) 
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and 



Thus the differential equation for H reduces to 

H + 2HH = H^H^ = AnGpH. (59) 

This coincides with the Hnear perturbation equation for 6f, (eq.[54]). Since H{t) 

is a decreasing function of t, this impHes that H(t) is the decaying sohition for 
equation (54). Then the corresponding growing solution D{t) can be obtained 
according to the standard procedure: subtracting equation (54) from equation 
(59) yields 



d 










m 


di 


dt 





0, (60) 



and therefore the formal expression for the growing solution in linear theory is 

It is often more useful to rewrite D{t) in terms of the redshift z as follows: 

Diz) = '-^Hiz)f^^dz', (62) 

where the proportional factor is chosen so as to reproduce D{z) 1/(1 + 2) for 
2; — > 00. Linear growth rates for the models described in §2.3 are summarized 
below. 

(a) Einstein - de Sitter model (f2„ = 1, = 0) 

Diz) = ^ (63) 

(b) Open model with vanishing cosmological constant (n„ < 1, VIa = 0) 



^W-l + ^3;i±^ln(./^T^-Vi), -^^^ (64) 
(c) Spatially-flat model with cosmological constant (0„ < 1,0^ = 1 — f^m) 
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For most purposes, the following fitting formulae [65] provide sufficiently accu- 
rate approximations. 



D{z) = 



1 



where 



2 VL^n{z) - \{z) + [1 + n{z)l2][l + A(z)/70] ' 

VL^^{\ + zf 



H{z) 



\{z) = 



Hp 
H{z) 



+ z)'^ + (1 - - r2A)(i + z)^ 
+ zf + (1 - a„ - f^A)(i + + ■ 



(66) 

(67) 



(68) 



(69) 



Note that ilm and J^a refer to the present values of the density parameter and 
the dimensionless cosmological constant, respectively, which will be frequently 
used in the rest of the review. 

Figure 2 shows the comparison of the numerically computed growth rate 
(thick lines) against the above fitting formulae (thin lines) which are practically 
indistinguishable . 
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3 Statistics of cosmological density fluctuations 



3.1 Gaussian random field 

Consider the density contrast Si = S{xi) = p{x)/p— 1 defined at the comoving 
position Xi- The density field is regarded as a stochastic variable, and thus forms 
a random field. The conventional assumption is that the primordial density field 
(in its linear regime) is Gaussian, i.e., its m-point joint probability distribution 

obeys the multi-variate Gaussian: 

P((5i, (52, ... , 5m) d6idd2 ■ ■ ■ ddm 



for an arbitrary positive integer to. Here Mij = (SiSj) is the covariancc matrix, 
and is its inverse. Since M^j = ^{xi,Xj), equation (70) implies that the 

statistical nature of the Gaussian density field is completely specified by the two- 
point correlation function ^ and its linear combination (including its derivative 
and integral). For an extensive discussion of the cosmological Gaussian density 
field, see ref. [5] . 

The Gaussian nature of the primordial density field is preserves in its linear 
evolution stage, but this is not the case in nonlinear stage. This is clear even 
from the definition of the Gaussian distribution: equation (70) formally assumes 
that the density contrast distributes symmetrically in the range of — oo < 5^ < 
DO, but in the real density field 6i cannot be less than —1. This assumption does 
not make any practical difference as long as the fluctuations are (infinitesimally) 
small, but is invalid in nonlinear regime where the typical amplitude of the 
fluctuations exceeds unity. 

In describing linear theory of cosmological density fluctuations, the Fourier 
transform of the spatial density contrast 6{x) = p{x)/{p) — 1 is the most basic 
variable: 




1 



exp 





(71) 




Then linear perturbation equation reads 




(73) 



(74) 
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Equation (74) yields <j){t) ex a '^{t)Dj^{t), and (j>{t) rapidly converges to a con- 
stant value. Thus Dj^ evolves following the growing solution in linear theory. 

The most popular statistic of clustering in the universe is the power spectrum 
of the density fluctuations: 

P{t,k)^{D^{tf) (75) 

which measures the amplitude of the mode of the wavcnumber fe. This is the 
Fourier transform of the two-point correlation function; 

^(a;, t) = -^ I P{t, fe)exp(-ife • x) dk. (76) 

Stt-^ J 

If the density field is globally homogeneous and isotropic (i.e., no preferred 
position or direction), equation (76) reduces to 

Since the above expression obtained after the ensemble average, x does not 
denote an amplitude of the position vector, but a comoving wavelength 27: /k 
corresponding to the wavenumber k = \k\. It should be noted that neither 
the power spectrum or the two-point correlation function contains information 
for the phase Thus in principle two clustering patterns may be completely 
different even if they have the identical two-point correlation functions. This im- 
plies the practical importance to describe the statistics of phases (/)^ in addition 
to the amplitude of clustering. 

In the Gaussian field, however, one can directly show that equation (70) 
reduces to the probability distribution function of and that are explicitly 
written as 

P{\5kl<t>kWkm - ||y-P (-^) d\5/-^. (78) 

mutually independently of k. The phase distribution is uniform, and thus does 
not carry information. The above probability distribution function is also de- 
rived when the real and imaginary parts of the Fourier components are 
uncorrelated and Gaussian distributed (with the dispersion P{k)/2) indepen- 
dently of k. As is expected, the distribution function (78) is completely fixed if 
P{k) is specified. This rephrases the previous statement that the Gaussian field 
is completely specified by the two-point correlation function in real space. 

Incidentally the one-point phase distribution turns out to be essentially uni- 
form even in a strongly non-Gaussian field [81, 21]. Thus it is unlikely to extract 
useful information directly out of it mainly due to the cyclic property of the 
phase. Very recently, however, Matsubara[52] and Hikage et al.[32] succeeded 
in detecting a signature of phase correlations in Fourier modes of mass density 
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fields induced by nonlinear gravitational clustering using the distribution func- 
tion of the phase sum of the Fourier modes for triangle wavevectors. Several 
different statistics which carry the phase information have been also proposed 
in cosmology, including the void probability function [95], the genus statistics 
[26], the Minkowski functional [58, 76]. 



3.2 Log- normal distribution 

A probability distribution function (PDF) of the cosmological density fluctua- 
tions is the most fundamental statistic characterizing the large-scale structure 
of the universe. As long as the density fluctuations are in the linear regime, 
their PDF remains Gaussian. Once they reach the nonlinear stage, however, 
their PDF significantly deviates from the initial Gaussian shape due to the 
strong non-linear mode-coupling and the non-locality of the gravitational dy- 
namics. The functional form for the resulting PDFs in nonlinear regimes are 
not known exactly, and a variety of phenomenological models have been pro- 
posed [34, 74, 9, 25]. 

Kayo et al. (2001) [40] showed that the one-point log-normal PDF: 



^S('5) = ^^exp 



{ln(l + ,5)+a?/2}2 



2at 



^ (79) 



1 + d 



describes very accurately the cosmological density distribution even in the non- 
linear regime (the rms variance (7„i ^ 4 and the over-density 6 ^ 100). The 
above function is characterized by a single parameter ai which is related to the 
variance of S. Since we use S to represent the density fluctuation field smoothed 
over R, its variance is computed from its power spectrum P^i explicitly as 

'^liiR) = TTT / Pnx{k)W^{kR)edk. (80) 

27r^ Jo 

Here we use subscripts "lin" and "nl" to distinguish the variables corresponding 
to the primordial (linear) and the evolved (nonlinear) density fields, respectively. 
Then a\ depends on the smoothing scale R alone and is given by 

al{R) = \n[l + al,{R)\- (81) 

Given a set of cosmological parameters, one can compute a'ni(i?) and thus ai{R) 
very accurately using a fitting formula for Pn\{k) (e.g., rcf. [65]). In this sense, 
the above log-normal PDF is completely specified without any free parameter. 

Figure 3 plots the one-point PDFs computed from cosmological N-body sim- 
ulations in SCDM, LCDM and OCDM (for Standard, Lambda and Open CDM) 
models, respectively [35, 40]. The simulations employ N = 256^ dark matter 
particles in a periodic comoving cube (100/i~^Mpc)^. The density fields are 
smoothed over Gaussian {Left panels) and Top-hat {Right panels) windows with 
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different smoothing lengths;i? — 2ft,^^Mpc, 6ft.^^Mpc and Solid lines 

show the log-normal PDFs adopting the value of cTni directly evaluated from sim- 
ulations (shown in each panel). The agreement between the log-normal model 
and the simulation results is quite impressive. A small deviation is noticeable 
only for (5 ^ — 0.5. 




1 00 



Figure 3: One-point PDFs in CDM models with Gaussian {left panels) and top- 
hat {right panels) smoothing windows; R = 2/i~^Mpc {cyan), 6/i~^Mpc {red), 
and 18/i~^Mpc {green). The solid and long- dashed lines represent the log-normal 
PDF adopting a^i calculated directly from the simulations and estimated from 
the nonlinear fitting formula of ref . [65], respectively, (ref. [40]) 



From an empirical point of view, Hubble (1934) [34] first noted that the 
galaxy distribution in angular cells on the celestial sphere may be approximated 
by a log-normal distribution, rather than a Gaussian. Theoretically the above 
log-normal function may be obtained from the one-to-one mapping between the 
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linear random-Gaussian and the nonlinear density fields [9]. Define a linear 
density field g smoothed over R obeying the Gaussian PDF: 



where the variance is computed from its linear power spectrum: 

1 

al^{R) = ^j^ P,,^{k)W\kR)edk. (83) 
If one introduces a new field 5 from g as 



1 + 5 = 




the PDF for 5 is simply given by {dg/dS)PQ (g) which reduces to equation (79). 

At this point, the transformation (84) is nothing but a mathematical pro- 
cedure to relate the Gaussian and the log-normal functions. Thus there is no 
physical reason to believe that the new field 5 should be regarded as a non- 
linear density field evolved from g even in an approximate sense. In fact it is 
physically unacceptable since the relation, if taken at face value, implies that 
the nonlinear density field is completely determined by its linear counterpart 
locally. We know, on the other hand, that the nonlinear gravitational evolution 
of cosmological density fluctuations proceeds in a quite nonlocal manner, and 
is sensitive to the surrounding mass distribution. Nevertheless the fact that 
the log-normal PDF provides a good fit to the simulation data empirically im- 
plies that the transformation (84) somehow captures an important aspect of the 
nonlinear evolution in the real universe. 



3.3 Higher-order correlation functions 

One of the most direct methods to evaluate the deviation from the Gaussianity is 
to compute the higher-order correlation functions. Suppose that Xi now labels 
the position of the i-th object (galaxy), and then the two-point correlation 
function ^12 = £^{xi,X2) is defined also in terms of the joint probability of the 
pair of objects located in the volume elements of SVi and 6V2. 

dPi2 = n^6Vi6V2 [1 + 62] . (85) 

where n is the mean number density of the objects. This definition is generalized 
to three- and four-point correlation functions, (123 = C(^i) ^2,353) and 771234 = 
77(3:1, X2,xs, X4), in a straightforward manner: 

SPi23= n^5Vi6V2SV3 [1 + 62 + 63 + 61 + Cm] , (86) 
SP1234 = n'^SVi5V25V36V4, [l + 62 + ^3 + ^4 + 63 + 64 + 64 
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+Cl23 + Cl24 + Cl34 + C234 

+6264 + 6364 + 6463 + ??1234] , (87) 

Apparently ^12, Ci23) and 771234 are symmetric with respect to the change of the 
indices. Define the following quantities with the same symmetry properties: 

■2^123 = ^12^23+^21^13+^23^31, (88) 

^1234 = ^12^23^34 + ^23^34^41 + 64^41^12 + ^13^32^24 + ^32^24^41 

+ ^24^41^13 + C12C24C43 + ^24^43^31 + ^31^12^24 + ^13^34^42 

+ 64^42^1 +?426lCl3, (89) 

-B1234 = ^2^13^4 + ^21^23^24 + ^31^32^34 + ^41^42^43, (90) 

C1234 = Cl23(Cl4 + C24 + C34) + Cl34(Cl2 + ^32 + ^42) 

+ Cl24(6l+62+64) + C234(a2+a3 + a4). (91) 

Then it is not unreasonable to suspect that the following relations hold 

C123 = Q ^123, (92) 
?7l234 = Ra ^1234 + -Rb -Bl234, (93) 

?7i234 = Rc C1234, (94) 

where Q. i?„. i?;,, and Rc are constants. In fact, the analysis of the two- 
dimensional galaxy catalogues [68] revealed 

Q = 1.29 ±0.21 (0.1/i-iMpc < r < 10/i-iMpc), (95) 
i?a = 2.5±0.6, Eb=4.3±1.2 (0.5/i-iMpc < r < 4/i-iMpc). (96) 

The generlization of those relations for A/'-point correlation functions is sus- 
pected to be hold generally: 

Af-l 

(ab) 

and called the hierarchical clustering ansatz. Cosmological N-body simulations 
approximately support the validity of the above ansatz, but also detect the finite 
deviation from it [82]. 

3.4 Genus statistics 

A complementary approach to characterize the clustering of the universe beyond 
the two-point correlation functions is the genus statistics [26] . This is a math- 
ematical measure of the topology of the isodensity surface. For definiteness, 
consider the density contrast field 6{x) at the position x in the survey volume 
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Vgji. This may be evaluated, for instance, by taking the ratio of the number of 
galaxies N{x, Vf) in the volume Vf centered at x to its average value N{Vf): 

Six, Vf) = ^^^^ - 1, aiVf) = {\S{x, Vf)\Y^^ (98) 
N(Vf) 

where (j{Vf) is its rms value. Consider the isodensity surface parameterized 
by a value of ly = 5{x,Vf)/(T{Vf). Genus is one of the topological numbers 
characterizing the surface defined as 



(99) 



where k is the Gaussian curvature of the isolated surface. The Gauss - Bonnet 
theorem implies that the value of g is indeed an integer and equal to the number 
of holes minus 1. This is qualitatively understood as follows; expand an arbitrary 
two-dimensional surface around a point as 

z = iox^ + bxy + ^cy^ = ^Kixj + ^K2xI. (100) 

Then the Gaussian curvature of the surface is defined hy k = kiK2- A surface 
topologically equivalent to a sphere (a torus) has k = 1 (k = 0), and thus 
equation (99) yields g = —1 {g — 0) which coincides with the number of holes 
minus 1. 

In reality, there arc many disconnected isodensity surfaces for a given u, and 
thus it is more convenient to define the genus density in the survey volume V^n 
using the additivity of the genus: 

0M^^,±9. = -J^^tl^^A. ,101) 

1 ^— 1 

where Ai {i = 1 ^ I) denote the disconnected isodensity surfaces with the same 
value oi V = 5{x,Vf)/(j{yf). Interestingly the Gaussian density field has an 
analytic expression for equation (101): 

G(.) = ^(^)'^'e-''^/^(l--^), (102) 

where 

(fc2) ^ ^ (103) 

/ P{k)W'^{kR)d^k 

is the moment of fc^ weighted over the power spectrum of fluctuations P{k), 
and the smoothing function W'^{kR) (e.g., ref. [5]). It should be noted that 
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in the Gaussian density field the information of the power spectrum shows up 
only in the proportional constant of equation (102), and its functional form is 
deterimined uniquely by the threshold value v. This z/-dependence reflects the 
the phase information which is ignored in the two-point correlation function and 
power spectrum. In this sense, genus statistics is a complementary measure of 
the clustering pattern of universe. 




Figure 4: Isodensity surfaces of dark matter distribution from N-body simula- 
tion; LCDM in (100/i-iMpc)3 sX v = -1.0, 0.0, 1.0, and 1.7. (ref. [54]) 

Even if the primordial density field obeys the Gaussian statistics, the subse- 
quent nonlinear gravitational evolution generates the significant non-Gaussianity. 
To distinguish the initial non-Gaussianity from that aquired by the nonlinear 
gravity is of fundamental importance in inferring the initial condition of the 
universe in a standard gravitational instability picture of structure formation. 
In a weakly nonlinear regime, Matsubara (1996) derived an analytic expression 
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for the non-Gaussianity emerging from the primordial Gaussian field [50]: 

I //J,2\\3/2 



where 



(104) 



(105) 



are the Hermite polynomials; Hi = i/, i?2 



i/^ — 61^^ + 3, i75 = i^'"' — lOi^^ + 15i^, The three quantities: 



((5^V^(5), t/ = -- 



;(V(5- V(5V2(5) (106) 



denote the third-order moments of ^. This expression plays a key role in un- 
derstanding if the non-Gaussianity in galaxy distribution is ascribed to the pri- 
mordial departure from the Gaussian statistics. 



3.5 Minkowski functionals 

In fact, gcmis is one of the complete sets of TV + 1 quantities, known as the 
Minkowski functionals (MFs), which determine the morphological properties of 
a pattern in AT-dimensional space. In the analysis of galaxy redshift survey data, 
one considers isodcnsity contours from the three-dimensional density contrast 
field (5 by taking its excursion set F^, i.e., the set of all points where the density 
contrast ^ exceeds the threshold level v as was the case in the case of genus 
described in the above subsection. 

All MFs can be expressed as integrals over the excursion set. While the first 
MF is simply given by the volume integration of a Heaviside step function 
normalized to the total volume Vtot , 

y^{y) = 7^ [ d'x@{iy - i^ix)) , (107) 
not Jv 

the other MFs, Vk{k = 1,2,3), are calculated by the surface integration of the 
local MFs, v^^^. The general expression is 

Vk{v) = ^l cPS{x)v'^%v,x), (108) 

I'tot JdF,, 

with the local Minkowski Functionals for k = 1, 2, 3 given by 

= i (109) 
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where i?i and R2 are the principal radii of curvature of the isodensity surface. 

For a 3-D Gaussian random field, the average MFs per unit volume can be 
expressed analytically as follows: 



where A = a = (S^)^/^, ai = (|V5p)i/2^ and 6 is the density 



The above MFs can be indeed interpreted as well-known geometric quanti- 
ties; the volume fraction Vb(z^), the total surface area Vi(i^), the integral mean 
curvature V2{i^), and the integral Gaussian curvature, i.e., the Euler character- 
istic V3{u). In our current definitions (eqs. [101] and [108], or eqs.[102] and 
[115]), one can easily show that V3{v) reduces simply to —G{h'). The MFs were 
first introduced to cosmological studies by Mecke et al. (1994) [58], and further 
details may be found in refs. [58, 31]. Analytic expressions of MFs in weakly 
non-Gaussian fields are derived in ref. [51]. 




(112) 



(113) 



(114) 



(115) 



contrast. 
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4 Galaxy biasing 



4.1 Concepts and definitions of biasing 

As discussed above, luminous objects, such as galaxies and quasars, are not 
direct tracers of the mass in the universe. In fact, the difference of the spatial 
distribution between luminous objects and dark matter, or the bias, has been 
indicated from a variety of observations. Galaxy "biasing" clearly exists. The 
fact that galaxies of different types cluster differently (e.g., ref. [16]) imphes that 
not all of them are exact tracers of the underlying mass distribution. See also 
§6- 

In order to confront theoretical model predictions for the mass distribution 
against observational data, one needs a relation of density fields of mass and 
luminous objects. The biasing of density peaks in a Gaussian random field is well 
formulated [37, 5] and it provides the first theoretical framework for the origin 
of galaxy density biasing. In this scheme, the galaxy-galaxy and mass-mass 
correlation functions are related in the linear regime via 

= b^mmir), (116) 

where the biasing parameter & is a constant independent of scale r. However, 
a much more specific linear biasing model is often assumed in common appli- 
cations, in which the local density fluctuation fields of galaxies and mass are 
assumed to be deterministically related via the relation 

Sg{X.) = bdmi^). (117) 

Note that equation (116) follows from equation (117), but the reverse is not 
true. 

The above deterministic linear biasing is not based on a reasonable physical 
motivation. If 6 > 1, it must break down in deep voids because values of 6g 
below —1 arc forbidden by definition. Even in the simple case of no evolution 
in comoving galaxy number density, the linear biasing relation is not preserved 
during the course of fluctuation growth. Non-linear biasing, where b varies with 
Sjn, is inevitable. 

Indeed, an analytical model for biasing of halos on the basis of the extended 
Press-Schechter approximation [59] predicts that the biasing in nonlinear and 
provides a useful approximation for its behavior as a function of scale, time and 
mass threshold. iV-body simulations provide a more accurate description of 
the nonlinearity of the halo biasing confirming the validity of the Mo & White 
model [36, 91]. 

4.2 Modeling biasing 

Biasing is likely to be stochastic, not deterministic [15]. An obvious part of this 
stochasticity can be attributed to the discrete sampling of the density field by 
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galaxies, i.e. the shot noise. In addition, a statistical, physical scatter in the 
efficiency of galaxy formation as a function of 6m is inevitable in any realistic 
scenario. For example, the random variations in the density on smaller scales is 
likely to be reflected in the efficiency of galaxy formation. As another example, 
the local geometry of the background structure, via the deformation tensor, 
must play a role too. Such 'hidden variables' would show up as physical scatter 
in the density-density relation [86]. 

Consider the density contrasts of visible objects and mass, 6„^j{x,z\R) and 
5„(x, z|i?), at a position x and a rcdshift z smoothed over a scale R [87]. In 
general, the former should depend on various other auxiliary variables A defined 
at different locations x' and redshifts z' smoothed over different scales R' in 
addition to the mass density contrast at the same position, 5„(x, While 
this relation can be schematically expressed as 

5obi(x, z\R) = ^[x, z, R, 5^(x, z\R),A{ic', z'\R'), ...], (118) 

it is impossible even to specify the list of the astrophysical variables A, and 
thus hopeless to predict the functional form in a rigorous manner. Therefore 
if one simply focuses on the relation between (5obj(x, ^;|i?) and 6^{x,z\R), the 
relation bcconicis incivitably stochastic and nonlinear due to the dependence on 
unspecified auxiliary variables A. 

For illustrative purposes, define the biasing factor as the ratio of the density 
contrasts of luminous objects and mass: 

B^,,{^,z\R)= ^_^(^^^|^) = 54x,.|i?) • (^^^^ 

Only in very idealized situations, the above nonlocal stochastic nonlinear factor 
in terms of may be approximated by 

1. a local stochastic nonlinear bias: 

B,,,(x, z\R) = 6<-.'[x, z, R, 5„(x, z\R), ^(x, z\R), . . .], (120) 

2. a local deterministic nonlinear bias: 

B,^,{^,z\R)=bltl'[z,R,S4^,z\R)], (121) 

and 

3. a local deterministic linear bias: 

S,bj(x,0|i?) = 6obj(^,i?) (122) 

From the above point of view, the local deterministic linear bias is obviously 
unrealistic, but is still a widely used conventional model for biasing. In fact, the 
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time- and scale-dependence of the linear bias factor 6obj(z, R) was neglected in 
many previous studies of biased galaxy formation until very recently. Currently, 
however, various models beyond the deterministic linear biasing have been se- 
riously considered with particular emphasis on the nonlinear and stochastic 
aspects of the biasing [71, 15, 86, 87]. 



4.3 Density peaks and dctrk matter halos as toy models 

for galcixy biasing 

Let us illustrate the biasing from numerical simulations by considering two 
specific and popular models; primordial density peaks and dark matter ha- 
los [87]. We use the N-body simulation data of _L = 100ft~"'^Mpc again for this 
purpose [35]. We select density peaks with the threshold of the peak height 
Vth = 1-0, 2.0, and 3.0. As for the dark matter halos, these are identified using 
the standard friend-of- friend algorithm with a linking length of 0.2 in units of 
the mean particle separation. We select halos of mass larger than the threshold 
Mth = 2.0 X 10^^, 5.0 X 10^2 and 1.0 x 10^^ h'^ Mq. 

Figures 5 and 6 depicts the distribution of dark matter particles (upper- 
panel), peaks {middle-panel) and halos (lower-panel) in LCDM model at z = 
and z = 2.2 within a circular slice (comoving radius of 150/i~^Mpc and thickness 
of 15h~^Mpc). We locate a fiducial observer in the center of the circle. Then 
the comoving position vector r for a particle with a comoving peculiar velocity 
V at a redshift z is observed at the position s in redshift space: 

1 r • V r 

s = r+-^^^, (123) 
H(z) |r| |r| 

where H(z) is the Hubble parameter at z. The right panels in Figures 5 and 

6 plot the observed distribution in redshift space, where the redshift-space dis- 
tortion is quite visible; the coherent velocity field enhances the structure per- 
pendicular to the line-of-sight of the observer (squashing) while the virialized 
clump becomes elongated along the line-of-sight (finger- of- God). 

We use two-point correlation functions to quantify stochasticity and nonlin- 
earity in biasing of peaks and halos, and explore the signature of the redshift- 
space distortion. Since we are interested in the relation of the biased objects and 
the dark matter, we introduce three different correlation functions; the auto- 
correlation functions of dark matter and the objects, ^,-„„ and and their 
cross-correlation function In the present case, the subscript o refers to ci- 
ther h (halos) or v (peaks) . Wc also use the superscripts R and S to distinguish 
quantities defined in real and redshift spaces, respectively. We estimate those 
correlation functions using the standard pair-coTint method. The correlation 
function ^'^^^ is evaluated under the distant-observer approximation. 

Those correlation functions are plotted in Figures 7 and 8 for peaks and 
halos, respectively. The correlation functions of biased objects generally have 
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larger amplitudes than those of mass. In nonlinear regimes > 1) the finger- 
of-god effect suppresses the amplitude of ^'^^ relative to while is larger 
than ^(^^ in linear regimes < 1) due to the coherent velocity field. 



(d) Real space Ked.sh.ift space z = o.O 




^ 300 h-iMpc 



Figure 5: Top- view of distribution of objects at z = in real (left panels) and 
redshift (right panels) spaces around the fiducial observer at the center; dark 
matter particles (top panels), peaks with ly > 2 (middle panels) and halos with 
M > 1.3 X IO^^Mq (bottom panels) in LCDM model. The thickness of those 
slices is IS/j-^Mpc. (ref. [87]) 
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Real space 



Kedshift space 



1 = 2.2 





Figure 7: Auto- and cross-correlation functions of dark matter and peaks in 
SCDM (left panels), LCDM (middle panels) and OCDM (right panels). Dif- 
ferent symbols indicate the results in real space (open squares for ly > 3, filled 
triangles for u > 2, open circles for v > 1, and crosses for dark matter), while 
different curves indicate those in redshift space (dashed for f > 3, dot-dashed 
for f > 2, solid for u > 1, and dotted for dark matter), (a) z = 0, (h) z = 2.2. 
(ref. [87]) 34 
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Figure 8: Same as Fig. 7 but for halo model; open squares and dashed lines for 
M > IO^^/j-^Mq, filled triangles and dot-dashed Hues for M > 5 x 10^2/i~^Mq, 
open circles and solid lines for M > 2 x lO^^/i^^M©, and crosses and dotted 
lines for dark matter. For SCDM model, we only plot the correlation functions 
with Mth = 5 X 10^2^ IO^^H-^Mq. (a) z = 0, (h) z = 2.2. (ref. [87]) 
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4.4 Biasing of galaxies in cosmological hydrodynamic sim- 
ulations 

Popular models of the biasing based on the peak or the dark halos are successful 
in capturing some essential features of biasing. None of the existing models 
of bias, however, seems to be sophisticated enough for the coming precision 
cosmology era. Development of a more detailed theoretical model of bias is 
needed. A straightforward next step is to resort to numerical simulations which 
take account of galaxy formation even if phenomenological at this point. We 
show an example of such approaches from Yoshikawa et al. (2001) [91] who apply 
cosmological smoothed particle hydrodynamic (SPH) )simulations in LCDM 
model with particular attention to the comparison of the biasing of dark halos 
and simulated galaxies (see also ref. [78]). 
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Figure 9: Distribution of gas particles, dark matter particles, galaxies and dark 
halos in the volume of 75/i~^ x 75h^^ x 30/i^^Mpc'^ model at z = 0. Upper- 
right:gas particles; Upper-left: dark matter particles; Lower-right: galaxies; 
Lower-left: DM cores, (ref. [91]) 
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Galaxies in their simulations are identified as clumps of cold and dense gas 
particles which satisfy the Jeans condition and have the SPH density more than 
100 times the mean baryon density at each redshift. Dark halos are identified 
with a standard friend-of-friend algorithm; the linking length is 0.164 times the 
mean separation of dark matter particles, for instance, at z — 0. In addition, 
they identify the surviving high-density substructures in dark halos, DM cores. 
See ref. [91] for further details. 

Figure 9 illustrates the distribution of dark matter particles, gas particles, 
dark halos and galaxies a,t z — where galaxies are more strongly clustered than 
dark halos. Figure 10 depicts a close-up snapshot of the most massive cluster at 
z = with mass M ~ 8 x 10^'' Mq. Circles in lower panels indicate the positions 
of galaxies identified in our simulation. 
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Figure 10: Snapshots of the most massive cluster {M ~ 8 x 10^** Mq) in the 
simulation at z = 0. Upper-left: dark matter; Upper-right: gas; Lower-left: DM 
cores; Lower-right: cold gas. Circles in lower panels indicate the positions of 
galaxies identified according to our criteria. The comoving size of the box is 
6.25/i-iMpc per side. (ref. [91]) 



Figure 11 shows the joint distribution of and 5g with mass density field 
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(5m at redshift z — 0, 1 and 2 smoothed over Rs = Uh^^Mpc. The conditional 
mean relation 6i{6m) computed directly from the simulation is plotted in solid 
lines while dashed lines indicate theoretical predictions of halo biasing by Taruya 
& Suto (2000) [86]. For a given smoothing scale, the simulated halos exhibit 
positive biasing for relatively small dm in agreement with the predictions. On the 
other hand, they tend to be underpopulated for large 6^^, or anti-biased. This is 
mainly due to the exclusion effect of dark halos due to their finite volume size 
which is not taken into account in the theoretical model. Since our simulated 
galaxies have smaller spatial extent than the halos, the exclusion effect is not so 
serious. This is clearly illustrated in lower panels in Figure 11, and indeed they 
show much better agreement with the theoretical model. 
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Figure 11: Joint probability distributions of overdensity fields for dark halos 
and galaxies with dark matter overdensity smoothed over Rs = 12/i~iMpc at 
redshift z = 0, 1 and 2. Solid lines indicate the conditional mean Si{6tn) for 
each object. Dashed lines in each panel depict the theoretical prediction of 
conditional mean by Taruya & Suto (2000). (ref. [91]) 



Turn next to a more conventional biasing parameter defined through the 
two-point statistics: 




(124) 



where ^ii{r) and ^mm(^) are two-point correlation functions of objects i and 
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of dark matter, respectively. While the above biasing parameter is ill-defined 
where either ^^(r) or £,mm{r) becomes negative, it is not the case at clustering 
scales of interest (< 10/i^^Mpc). 
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Figure 12: Two-point correlation functions of dark matter, galaxies, and dark 
halos from cosmological hydrodynamical simulations, (ref. [91]) 

Figure 12 shows two-point correlation functions of dark matter, galaxies, 
dark halos and DM cores {upper and middle panels), and the profiles of biasing 
parameters b^{r) for those objects {lower panels) at z = 0, 1 and 2. In lower 
panels, we also plot the parameter 6var,i = fi/cm, which are defined in terms 
of the one-point statistics (variance), for comparison on the smoothing scale 

= Ah^^Mpc, 8/i~^Mpc and 12ft,^^Mpc at r = Rs for each kind of objects 
by different symbols. In the upper panels, we show the correlation functions of 
DM cores identified with two different maximum linking length; Zmax — 0.05 and 
'max = &h/2. Correlation functions of DM cores identified with ^max = 0.05 are 
similar to those of galaxies. On the other hand, those identified with Zmax — &h/2 
exhibit much weaker correlation, and are rather similar to those of dark halos. 
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This is due to the fact that the present algorithm of group identification with 
larger /max tends to pick up lower mass halos which are poorly resolved in our 
numerical resolution. 

The correlation functions of galaxies are almost unchanged with redshift, 
and correlation functions of dark halos only slightly evolve between z = and 
2. By contrast, the amplitude of the dark matter correlation functions evolves 
rapidly by factor of ~ 10 from z = 2 to z = 0. The biasing parameter 6^_g 
is larger at a higher redshift, for example, b^^g ~ 2-2.5 at z = 2. The biasing 
parameter b^^h for dark halos is systematically lower than that of galaxies and 
DM cores again due to the volume exclusion effect. At z = 0, galaxies and DM 
cores are slightly anti-biased relative to dark matter at r ~ l/i^^Mpc. In lower 
panels, we also plot the one-point biasing parameter 6var,i = oija^ ai r = Rg 
for comparison. In general we find that b^A is very close to 5var,j at z ~ 0, but 
systematically lower than 6var,i at higher redshifts. 




Figure 13: Two-point correlation functions for the old and young populations 
of galaxies at z = as well as that of dark matter distribution. The profiles 
of bias parameters b^ (r) for both of the two populations are also shown in the 
lower panel, (ref. [91]) 



For each galaxy identified at z = 0, we define its formation redshift Zf by the 
epoch when half of its cooled gas particles satisfy our criteria of galaxy formation. 
Roughly speaking, Zf corresponds to the median formation redshift of stars in 
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the present-day galaxies. We divide all simulated galaxies at z = into two 
populations (the young population with Z{ < 1.7 and the old population with 
Z{ > 1.7) so as to approximate the observed number ratio of 3/1 for late- type 
and early- type galaxies. 

The difference of the clustering amplitude can be also quantified by their two- 
point correlation functions at z = as plotted in Figure 13. The old population 
indeed clusters more strongly than the mass, and the young population is anti- 
biased. The relative bias between the two populations 6^^g = -\/Coid/Cyoung 
ranges 1.5 and 2 for lft,~'^Mpc < r < 20/i~"'^Mpc, where ^young and ^oid are the 
two-point correlation functions of the young and old populations. 

4.5 Halo occupation function approach for galaixy biasing 

Since the clustering of dark matter halos is well understood now, one can de- 
scribe the galaxy biasing if the halo model is combined with the relation between 
the halos and luminous objects. This is another approach to galaxy biasing, halo 
occupation function (HOF), which becomes very popular recently. Indeed the 
basic idea behind HOF has a long history, but the model predictions have been 
significantly improved with the recent accurate models for the mass function, 
the biasing and the density profile of dark matter halos. We refer the readers 
for an extensive review on the HOF by Cooray & Sheth (2002) [12]. Here we 
briefly outline this approach. 

We adopt a simple parametric form for the average number of a given galaxy 
population as a function of the hosting halo mass: 



The above statistical and empirical relation is the essential ingredient in the 
current modeling characterized by the minimum mass of halos which host the 
population of galaxies (Mmin), a normalization parameter which can be inter- 
preted as the critical mass above which halos typically host more than one 
galaxy (Mi ; note that Mi may exceed Mmin since the above relation represents 
the statistical expected value of number of galaxies), and the power-law index 
of the mass dependence of the efficiency of galaxy formation (a). We will put 
constraints on the three parameters from the observed number density and clus- 
tering amplitude for each galaxy population. In short, the number density of 
galaxies is most sensitive to Mi which changes the average number of galaxies 
per halo. The clustering amplitude on large scales is determined by the hosting 
halos and thus very sensitive to the mass of those halos, Mmin- The cluster- 
ing on smaller scales, on the other hand, depends on those three parameters 
in a fairly complicated fashion; roughly speaking, Mmin changes the amplitude, 
while a, and to a lesser extent Mi as well, changes the slope. 

With the above relation, the number density of the corresponding galaxy 




(M > Mmin) 
(M<Mmin) ■ 



(125) 
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population at redshift z is given by 

f-OC 

ng,,{z)= dM n^^x^{M,z) N^{M), (126) 

■/M„i„ 

where nhaio(^'^) denotes the halo mass function. 

The galaxy two-point correlation function on small scales is dominated by 
contributions of galaxy pairs located in the same halo. For instance. Bullock 
et al. (2002) [8] adopted the mean number of galaxy pairs {Ng{Ng — 1))(M) 
within a halo of mass M of the form: 

r NliM) (7Vg(M) > 1) 

(iVg(iVg - 1))(M) = { iV|(M) log(4iVg(M))/log(4) (1 > iVg(Af) > 0.25) . 

[ (N^iM) < 0.25) 

(127) 

In the framework of the halo model, the galaxy power spectrum consists of 
two contributions, one from galaxy pairs located in the same halo (1-halo term) 
and the other from galaxy pairs located in two different halos (2-halo term): 

P|°*(fc) = Pgi''(A;) + Pf (fc). (128) 

The 1-halo term is written as 

= /9 L2 I dMn^^M) {Ng{N,-l))b{MMk,M)\P. (129) 
yzTV) Tig 2 J 

Seljak (2000) [77] chose p = 2 for {Ng{Ng - 1)) > 1 and p = 1 for (iVg(A^g - 1)) < 
1. The 2-halo term on the assumption of the linear halo bias model [59] reduces 
to 



Plin(fe) 



J dM n^^xo{M)N^{M)b{M)y{k,M) 



(130) 



where Piin(fc) is the linear dark matter power spectrum, h{M) is the halo bias 
factor, and M) is the Fourier transform of the halo dark matter profile 
normalized by its mass, y{k,M) — p{k,M)/M [77]. 

The halo occupation formalism, although simple, provides a \iseful frame- 
work in deriving constraints on galax;y formation models from large data sets of 
the upcoming galaxy redshift surveys. For example, Zehavi et al. (2003) [99] 
used the halo occupation formalism to model departures from a power law in the 
SDSS galaxy correlation function. They demonstrated that this is due to the 
transition from a large-scale regime dominated by galaxy pairs in different halos 
to a small-scale regime dominated by those in the same halo. Magliocchetti & 
Porciani (2003) [47] applied the halo occupation formalism to the 2dFGRS clus- 
tering results per spectral type of Madgwick et al. (2003) [45]. This provides 
constraints on the distribution of late-type and early-type galaxies within the 
dark matter halos of different mass. 
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5 Relativist ic effects observable in clustering at 
high redshifts 

Redshift surveys of galaxies definitely serve as the central database for obser- 
vational cosmology. In addition to the existing shallower surveys {z < 0.2), 
clustering in the; universe in the range z = 1 — 3 has been partially revealed by, 
for instance, the Lyman-break galaxies and X-ray selected AGNs. In particu- 
lar, the 2dF and SDSS QSO redshift surveys promise to extend the observable 
scale of the universe by an order of magnitude, up to a few Gpc. A proper 
interpretation of such redshift surveys in terms of the clustering evolution, how- 
ever, requires an understanding of many cosmological effects which can be ne- 
glected for z <C 1 and thus have not been considered seriously so far. These 
cosmological contaminations include linear redshift-space (velocity) distortion, 
nonlinear redshift-space (velocity) distortion, cosmological redshift-space (geo- 
metrical) distortion, and cosmological light-cone effect. 

We describe a theoretical formalism to incorporate those eff'ects, in particular 
the cosmological redshift-distortion and light-cone effects, and present several 
specific predictions in CDM models. The details of the material presented in 
this section may be found in refs. [83, 89, 90, 46, 28, 29]. 

5.1 Cosmological lightcone effect on the two-point corre- 
lation functions 

Observing a distant patc;h of the universe is equivalent to observing the past. 
Due to the finite light velocity, a line-of-sight direction of a redshift survey is 
along the time, as well as spatial, coordinate axis. Therefore the entire sample 
does not consist of objects on a constant-time hypersurface, but rather on a light- 
cone, i.e., a null hypersurface defined by observers at 2 = 0. This implies that 
many properties of the objects change across the depth of the survey volume, 
including the mean density, the amplitude of spatial clustering of dark matter, 
the bias of luminous objects with respect to mass, and the intrinsic evolution of 
the absolute magnitude and spectral energy distribution. These aspects should 
be properly taken into account in order to extract cosmological information from 
observed samples of redshift surveys. 

In order to predict quantitatively the two-point statistics of objects on the 
light cone, one must take account of (i) nonlinear gravitational evolution, (ii) 
linear redshift-space distortion, (iii) nonlinear redshift-space distortion, (iv) 
weighted averaging over the light-cone, (v) cosmological redshift-space distor- 
tion due to the geometry of the universe, and (vi) object-dependent clustering 
bias. The efi'ect (v) comes from our ignorance of the correct cosmological pa- 
rameters, and (vi) is rather sensitive to the objects which one has in mind. Thus 
the latter two effects will be disciisscd in the next subsection. 

Nonlinear gravitational evolution of mass density fluctuations is now well 
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understood, at least for two-point statistics. In practice, we adopt an accurate 
fitting formula [65] for the nonlinear power spectrum P^{k, z) in terms of its 
linear counterpart. If one assumes a scale-independent deterministic linear bias, 
furthermore, the power spectrum distorted by the peculiar velocity field is known 
to be well approximated by the following expression: 



P^^\k^,ku;z) = b''iz)Pi^l{k;z) 



L»vei [fc||ap(^)] , (131) 



where k± and k\\ are the comoving wavenumber perpendicular and parallel to the 
line-of-sight of an observer, and Pj^}^{k; z) is the mass power spectrum in real 
space. The second factor in the right-hand-side comes from the linear redshift- 
space distortion [38], and the last factor is a phenomenological correction for 
non-linear velocity effect [65]. In the above, we introduce 
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We assume that the pair-wise velocity distribution in real space is approxi- 
mated by 

1 / V2\V^2\\ 



fv{vi2) 



V2, 



exp 



(Tp 



(133) 



with 0-p being the 1-dimensional pair-wise peculiar velocity dispersion. Then 
the finger-of-god effect is modeled by the damping function, Dvei [A;||(Tp(2)] : 



Dvcl[fc/iO-p 



1 



1+k2^2^ 



(134) 



where /i is the direction cosine in fc-spacc, and the dimensionless wavemimber 
K is related to the peculiar velocity dispersion a-p in the physical velocity units: 



k{z) = 



k{l+z)ap{z) 
V2H{z) 



(135) 



Since we are mainly interested in the scales around l/i~^Mpc, we adopt the 
following fitting formula throughout the analysis below which better approxi- 
mates the small-scale dispersions in physical units: 



(Tp{z) 



740(1 + z)-ikm/s 
650(1 + z)-"-*km/s 



for SCDM model 
for ACDM model. 



(136) 



Integrating equation (131) over fj,, one obtains the direction-averaged power 
spectrum in redshift space: 



A{K) + ^(3{z)B{i^) + ^(3\z)C{n) 



(137) 
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where 



A{k) = 
Bin) = 



C{k) = 



arctan(K) 



1 - 



arctan(K) 



3k2 



3 3 arctan(K) 

1 o ~t~ o 



(138) 
(139) 
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Adopting those approximations, the direction-averaged correlation functions 
on the hght-cone are finally computed as 



[(t){z)no{z)fS,{xs;z) 



[(l){z)no{z)f 



(141) 



where Zmin and Zmax denote the redshift range of the survey, and 



(142) 



Throughout the present analysis, we assume a standard Robertson - Walker 
metric of the form: 



= _dt^ + a{tf{dx^ + SKixflde"^ + sin^ Odcj)'^]), 
where Sk (x) is determined by the sign of the curvature K as 



sin {y/Kx)HK {K > 0) 

Sk{x) = { X {K = Q) 

sinh (\/^x) /\/^ {K < 0) 



(143) 



(144) 



where the present scale factor Oo is normalized as unity, and the spatial curvature 
K is given as (eq.[13]) 

K = Hl{n^ + n^~i). (145) 

The radial comoving distance x{z) is computed by 



X{z 



dt 



, a{t) 7o ii{z) 



(146) 



The comoving angular diameter distance Dciz) at redshift z is equivalent to 
5~^(x(^)), and, in the case of Oa = 0, is explicitly given by Mattig's formula: 



D,{z) = 



l + z + y/T+JTZz 



Hol + zl + n^z/2 + ^/l + Q„z' 



(147) 
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Then dVc/dz, the comoving volume element per unit solid angle, is explicitly 
given as 

= ^^^^^ (148) 

Ho^n^il + z)3 + (1 - n„ - + zf + 

5.2 Evaluating two-point correlation functions from N- 
body simulation data 

The theoretical modeling described above was tested against simulation results 
by Hamana, Colombi & Suto (2001) [28]. Using cosmological A''-body simula- 
tions in SCDM and ACDM models, they generated light-cone samples as follows; 
first, they adopt a distance observer approximation and assume that the line-of- 
sight direction is parallel to Z-axis regardless with its (X, F) position. Second, 
they periodically duplicate the simulation box along the Z-dircction so that 
at a redshift 2, the position and velocity of those particles locating within an 
interval x(-2) ± ^xi.^^ a^re dumped, where Ax(-2) is determined by the output 
time-interval of the original A^-body simulation. Finally they extract five in- 
dependent (non-overlapping) cone-shape samples with the angular radius of 1 
degree (the field-of-view of tt degree^). In this manner, they have generated 
mock data samples on the light-cone continuously extending up to z = 0.4 (rel- 
evant for galaxy samples) and z = 2.0 (relevant for QSO samples), respectively 
from the small and large boxes. 

Two-point correlation function is estimated by the conventional pair-count 
adopting the following estimator [43]: 

c(.) = (149) 

ixix\x) 

The comoving separation xn of two objects located at z\ and zi with an 
angular separation 6\i is given by 

^12 = x\+ xi^ — Kx\x\{\ -\- co^ Q\'2) 



-2xiX2\Jl- Kxl-^l- Kxlcosei2, (150) 

where xi = Dc{zi) and X2 = Dc{z2)- 

In redshift space, the observed redshift ^;obs for each object differs from the 
"real" one ^reai due to the velocity distortion effect: 

•2obs — -^real "I" (1 + ■2'real)^pec) (-'-^-'-) 

where t^pec is the line of sight relative peculiar velocity between the object and 
the observer in physical units. Then the comoving separation si2 of two objects 
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redshift space is computed as 

3l + S2 - KsfsKl + COS^ 012) 



.2 

*12 



= si 



— 2siS2\Jl — Ks\^ 1 — cos 6*12, 

where si = i?c(2obs,i) and S2 = £'c(2obs,2)- 
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Figure 14: Mass two-point correlation functions on the light cone for particles 
with redshift-dependent selection functions in SCDM model. Upper: z < 0.4, 
Lower: 0.2 < z < 2.0. Left: with selection function whose shape is the same as 
that of the B-band magnitude limit of 19 for galaxies (upper) and 21 for QSOs 
(lower). Right: randomly selected N ~ 10^ particles from the particles in the 
left results, (ref. [28]) 



In properly predicting the power spectra on the light cone, the selection func- 
tion should be specified. For galaxies, they adopt a B-band luminosity function 
of the APM galaxies fitted to the Schechter function [44]. For quasars, they 
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adopt the B-band luminosity function from the 2dF QSO survey data [7]. To 
compute the B-band apparent magnitude from a quasar of absolute magnitude 
Mb at z (with the luminosity distance d\,{zy)^ they applied the K-correction: 

B = A/b + 5 log((iL(z)/10pc) - 2.5(1 - v) log(l + z) (153) 

for the quasar energy spectrum L^, cc v~'p (they used p = 0.5). In practice, 
they adopt the galaxy selection function (/)gai(< -Biim,^) with — 19 and 
2^inin = 0.01 for the small box realizations, and the QSO selection function 
0QSo(< B\{a\,z) with Slim — 21 and z,nin — 0.2 for the large box realizations. 
They do not introduce the spatial biasing between selected particles and the 
underlying dark matter. 

Figures 14 and 15 plot the two-point correlation functions in SCDM and 
ACDM, respectively, taking account of the selection functions. It is clear that 
the simulation results and the predictions are in good agreement. 




Figure 15: Same as Figure 14 but for ACDM model, (ref. [28]) 
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5.3 Cosmological redshift-space distortion 

Consider a spherical object at high redshift. If the wrong cosmology is as- 
sumed in interpreting the distance-redshift relation along the line of sight and 
in the transverse direction, the sphere will appear distorted. Alcock & Paczyn- 
ski (1979) [2] pointed out that this curvature effect could be used to estimate 
the cosmological constant. Matsubara & Suto (1996) [54] and Ballinger, Pea- 
cock & Heavens (1996) [4] developed a theoretical framework to describe the 
geometrical distortion effect (cosmological redshift distortion) in the two-point 
correlation function and the power spectrum of distant objects, respectively. 
Certain studies were less optimistic than others about the possibility of measur- 
ing this AP effect. For example, Ballinger, Peacock and Heavens (1996) argued 
that the geometrical distortion could be confused with the dynamical redshift 
distortions caused by peculiar velocities and characterized by the linear theory 

parameter (3 = — g^. Matsubara & Szalay (2002, 2003) [55, 56] showed that the 
typical SDSS and 2dF samples of normal galaxies at low redshift (~ 0.1) have 
sufficiently low signal-to-noise, but they are too shallow to detect the Alcock & 
Paczynski effect. On the other hand, the quasar SDSS and 2dFGRS surveys are 
at a useful redshift, but they are too sparse. A more promising sample is the 
SDSS Luminous Red Galaxies survey (out to redshift z ^ 0.5) which turns out 
to be optimal in terms of both depth and density. 

While this analysis is promising, it remains to be tested if non-linear cluster- 
ing and complicated biasing (which is quite plausible for red galaxies) would not 
'contaminate' the measurement of the Equation of State. Even if the Alcock & 
Paczynski test turns out to be less accurate than other cosmological tests (e.g., 
CMB and SN la) the effect itself is an interesting and important ingredient in 
analyzing the clustering pattern of galaxies at high redshifts. We shall now 
present the formalism for this effect. 

Due to a general-rclativistic effect through the geometry of the universe, the 
observable separations perpendicular and parallel to the line-of-sight direction, 
a;sj. = {c/Ho)zd6 and Xs\\ = {c/Ho)Sz, are mapped differently to the correspond- 
ing comoving separations in real space x± and x^: 

Xs±{z) = x±cz/[Ho{l + z)dj^{z)]=x±/c±{z), (154) 
Xsn{z) = xiiH{z)/Ho = xii/cii{z), (155) 

with dA{z) being the angular diameter distance. The difference between c±{z) 
and C|| (z) generates an apparent anisotropy in the clustering statistics, which 
should be isotropic in the comoving space. Then the power spectrum in cosmo- 
logical redshift space, 

p(cRD)^ is related to P(^) defined in the comoving redshift 

space as 

P^'^--\h^,h\i;z) = _-L^p(S) ,) , (156) 

C±{zycii{z) \C±{Z) C||(^) J 
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where the first factor comes from the Jacobian of the volume element dk'^^dksw, 
and ksi_ = c±{z)k± and fcs|| — C||(z)fc|| are the wavenumber perpendicular and 
parallel to the line-of-sight direction. 
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Figure 16: Two-dimensional power spectra in cosmological redshift space at 
z = 2.2. (ref. [46]) 



Using equation (131), equation (156) reduces to 
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Figure 16 shows anisotropic power spectra P'^'^^°\ks, ^ik] z — 2.2). As spe- 
cific examples, we consider SCDM, LCDM and OCDM models, which have 
Oa, /i, crs) = (1.0, 0.0, 0.5, 0.6), (0.3, 0.7, 0.7, 1.0), and (0.3, 0.0, 0.7, 1.0), re- 
spectively. Clearly the linear theory predictions (cTp = 0; top panels) are quite 
different from the results of N-body simulations (bottom panels) , indicating the 
importance of the nonlinear velocity effects (cp computed according to ref . [57] ; 
middle panels). 




0.5 0.5 0.5 1 



% 

Figure 17: The confidence contours on ri„,-ilA plane from the x^-analysis of the 
monopole and quadrupole moments of the power spectrum in the cosmological 
redshift space at z = 2.2. We randomly selected TV = 5 x 10^ (upper panels), 
= 5 X 10"* (middle panels), and = 5 x 10^ (lower panels) particles from 
N-body simulation. The value of cts is adopted from the cluster abundance, 
(ref. [46]) 
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Next we decompose the power spectrum into harmonics: 



P{k,Hk\z) = Pi{k)Li{tik), Pi{k;z) 




where Li[fj,k) are the l-ih order Legendre polynomials. Similarly, the two-point 
correlation function is decomposed as 



using the direction cosine, fix, between the separation vector and the line-of- 
sight. The above multipole moments satisfy the following relations: 



with ji{kx) being spherical Bessel functions. Substituting P('^^°)(A;s, /x^; z) in 
equation (159) yields P[''''°\h; z) , and then^(°^°)(xs;z) can be computed from 
equation (161). 

Comparison of the monopoles and quadrupoles from simulations and model 
predictions exhibits how the results are sensitive to the cosmological parame- 
ters, which in turn may put potentially useful constraints on (ri^,^]^). Figure 
17 indicates the feasibility, which interestingly results in a constraint fairly or- 
thogonal to that from the Supernovae la Hubble diagram. 

5.4 Two-Point Clustering Statistics on a Light-Cone in 
Cosmological Redshift Space 

In order to explore the relation between the two-point statistics on a constant- 
time hypersurface in real space and that on a light-cone hypersurface in cosmo- 
logical redshift space, we simply consider the case of the deterministic, linear, 
and scale- independent bias: 



In what follows, we explicitly use the subscript, mass, to indicate the quantities 
related to the mass density field, and those without the subscript correspond to 
objects satisfying equation (163). 

Using equation (157), the two-point correlation function in the cosmological 
redshift space, ^^^^°\xs±,Xsi\; z), is computed as 





(162) 



(161) 



(5(x, z) = b{z) 6yn{x,z). 



(163) 
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= ^yp(s)(k;^)expHk.x)d3fc 

= ^(^)(c_La;sj_,C||a;s||;2;), (164) 

where ^'^'(x^. .T||; 2;) is the redshift-space correlation function defined through 
equation (131). 

bmce Pt'"'\k,;z) andej""°'(x.;z) are defined in redshift space, the proper 
weight should be 

^3^(cRD) [0(^)^cRD^^^j2 ^ d3.T[<^(z)no™(z)]2 cx(2)2c||(z), (165) 

where nQ^°{z) and nQ°'"(z) denote the number densities of the objects in cosmo- 
logical redshift space and comoving space, respectively, and (/){z) is the selection 
function determined by the observational target selection and the luminosity 
function of the objects. Then, the final expressions [84] reduce to 

dz^ [cf>{z)nr{z)]W{z)\\{z)P[''''''\h;z) 

/ ^^~d^ [^^(^)"o°""(2)]^c_l(2)2c||(2;) 

/ dz-^ [^{z)nr{z)]W{zfc^^{z)^r-{x,;z) 

^^°'°'°'(^s) = JJ ,(167) 

/ '^^H^ y>{z)nQ"'{z)fc^{zfc\\{z) 

where ^min and 2;max denote the redshift range of the survey, dVc /dz = d"^ (z) /H{z) 
is the comoving volume element per unit solid angle. 

Note that fcg and Xg, defined in P/°^°^(fcs;z) and ^"(cCg; z), are related 
to their comoving counterparts at z through equations (158) and (154), while 
those in p/^'^'*^^°^(fcs) and ^i^'^'^^°\xs) are not specifically related to any co- 
moving wavenumber and separation. Rather, they correspond to the quan- 
tities averaged over the range of z satisfying the observable conditions Xg = 
{c/Ho)VSz^ + z'^Se'^ and kg = 27r/xg. 

Let us show specific examples of the two-point clustering statistics on a light- 
cone in cosmological redshift space. We consider SCDM and LCDM models, 
and take into account the selection functions relevant to the upcoming SDSS 
spectroscopic samples of galaxies and quasars by adopting the B-band limiting 
magnitudes of 19 and 20, respectively. 

Figure 18 compares the predictions for the angle-averaged (monopole) power 
spectra under various approximations. The upper and lower panels adopt the 
selection functions appropriate for galaxies in < 2; < Zmax = 0.2 and QSOs in 
< 2; < 2;max = 5, respectively. The left and right panels present the results in 
SCDM and LCDM models. For simplicity we adopt a scale-independent linear 
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bias model [23]: 

6(z) = l + -^[6(fc,z = 0)-l], (168) 

with b{k, z = 0) = 1 and 1.5 for galaxies and quasars, respeetively. 

The upper and lower panels correspond to magnitude-limited samples of 
galaxies < 19 in < z < Zmax — 0.2; no bias model) and QSOs {B < 20 in 
< z < Zmax = 5; Fry's linear bias model), respectively. We present the results 
normalized by the real-space power spectrum in linear theory z) [5], 

and P^^\k;z = 0), P^^\k;z = Zmax), PQ'^'^°\ks; z = Zmax), and P^'''^''^'^°\ks) 
are computed using the nonlinear power spectrum [65]. 
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Figure 18: Light-cone and cosmological redshift-space distortion effects on 
angle-averaged power spectra, (ref. [84]) 



Consider first the results for the galaxy sample (upper panels). On linear 
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scales {k < O.lh Mpc~-^), Pq {k;z = 0) plotted in dashed lines is enhanced rel- 
ative to that in real space, mainly due to a linear redshift-space distortion [the 
Kaiser factor in equation (131)]. For nonlinear scales, the nonlinear gravitational 
evolution increases the power spectrum in real space, while the finger-of-god ef- 
fect suppresses that in redshift space. Thus, the net result is sensitive to the 
shape and the amplitude of the fluctuation spectrum, itself; in the LCDM model 
that we adopted, the nonlinear gravitational growth in real space is stronger 
than the suppression due to the finger-of-god effect. Thus, Pq {k;z = 0) be- 
comes larger than its real-space counterpart in linear theory. In the SCDM 
model, however, this is opposite and Pq {k;z = 0) becomes smaller. 

The power spectra at 2; = 0.2 (dash-dotted lines) are smaller than those at 
2: = by the corresponding growth factor of the fluctuations, and one might ex- 
pect that the amplitude of the power spectra on the light-cone (solid lines) 
would be in-between the two. While this is correct, if we use the comov- 
ing wavcnumbcr, the actual observation on the light-cone in the cosmologi- 
cal redshift space should be expressed in terms of kg [equation (158)]. If we 
plot the power spectra a,t z = 0.2 taking into account the geometrical distor- 
tion, PQ^^°\ks;z = 0.2) in the dotted lines becomes significantly larger than 
PQ^\k; z = 0.2). Therefore, Plf^^'^^°\ks) should take a value between those of 
Pt'^°\kg;z = 0) = P^^\k;z = 0) and P^'''^°\kg; z = 0.2). This explains the 
qualitative features shown in the upper panels of Figure 18. As a result, both 
the cosmological redshift-space distortion and the light-cone effect substantially 
change the predicted shape and amplitude of the power spectra, even for the 
galaxy sample [60] . The results for the QSO sample can be basically understood 
in a similar manner, except that the evolution of the bias makes a significant 
difference, since the sample extends to much higher redshifts. 

Figure 19 shows the results for the angle-averaged (monopole) two-point 
correlation fimctions, exactly corresponding to those in Figure 18. The results 
in this figure can also be understood by an analogy of those presented in Figure 
18 at A; ~ 27r/a;. Unlike the power spectra, however, two-point correlation 
functions are not positive deflnite. The funny features in Figure 19 on scales 
larger than SO/j-^Mpc (lOOft-^Mpc) in SCDM (LCDM) originate from the fact 
that ^^^'^™\x,z = 0) becomes negative there. 

In fact, since the resulting predictions are sensitive to the bias, which is 
unlikely to quantitatively be specified by theory, the present methodology will 
find two completely different applications. For relatively shallower catalogues, 
like galaxy samples, the evolution of bias is not supposed to be so strong. Thus, 
one may estimate the cosmological parameters from the observed degree of the 
redshift distortion, as has been conducted conventionally. Most importantly, we 
can correct for the systcmatics due to the light-cone and geometrical distortion 
effects, which affect the estimate of the parameters by ~ 10%. Alternatively, for 
deeper catalogues like high-redshift quasar samples, one can extract information 
on the object-dependent bias only by correcting the observed data on the basis 
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Figure 19: Same as Figure 18 on angle-averaged two-point correlation functions, 
(ref. [84]) 



of our formulae. 

In a sense, the former approach uses the light-cone and geometrical dis- 
tortion effects as real cosmological signals, while the latter regards them as 
inevitable, but physically removable, noise. In both cases, the present method- 
ology is important in properly interpreting the observations of the universe at 
high redshifts. 
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6 Recent results from 2dF and SDSS 



6.1 The latest galaxy redshift surveys 

Redshifts surveys in the 1980s and the 1990s (e.g the CfA, IRAS and Las cam- 
panas surveys) measured thousands to tens of thousands galaxy redshifts. Mul- 
tifibre technology now allows us to measure redshifts of millions of galaxies. 
Below we summarize briefly the properties of the main new surveys, 2dFGRS, 
SDSS, 6dF, VIRMOS, DEEP2 and we discuss key results from 2dFGRS and 
SDSS. Further analysis of these surveys is currently underway. 

6.1.1 The 2dF galaxy redshift survey 

The Anglo- Australian 2 degree Field Galaxy Redshift (2dFGRS)^ has recently 
been completed with redshifts for 230,000 galaxies selected from the APM cat- 
alogue December 2002) down to an extinction corrected magnitude limit of 
bj < 19.45. The main survey regions are two declination strips, in the northern 
and southern Galactic hemispheres, and also 100 random fields, covering in to- 
tal about 1800 deg2 (Figs. 20 and 21). The median redshift of the 2dFGRS is 
z ~ 0.1 (see refs. [11, 67] for reviews). 

Harm Pdg 




_*PH mnrwi FiH* 

Figure 20: The 2dFGRS fields (small circles) superimposed on the APM cata- 
logue area (dotted outlines of Sky Survey plates). 



http; / / www.mso. anu.edu. au / 2dFGRS / 
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Figure 21: The distribution of 63,000 2dFGRS galaxies in the NGP (Left) and 
SGP {Right) strips. 



6.1.2 The SDSS galaxy redshift survey 

The SDSS (Sloan Digital Sky Survey) is a US- Japan-Germany joint project to 
image a quarter of the Celestial Sphere at high Galactic latitude as well as to 
obtain spectra of galaxies and quasars from the imaging data^. The dedicated 
2.5 meter telescope at Apache Point Observatory is equipped with a multi-CCD 
camera with five broad bands centered at 3561, 4676, 6176, 7494, and 8873 A. 
For further details of SDSS, see refs. [97, 80] 

The latest map of the SDSS galaxy distribution, together with a typical slice, 
are shown in Figures 22, and 23 [31]. The three-dimensional map centered on 
us in equatorial coordinate system is shown Figure 22. Redshift slices of galaxies 
centered around the equatorial plane with various redshift limits and thicknesses 
of planes are shown in Figure 23: z < 0.05 with thickness of 10/i~^Mpc centered 
around the equatorial plane in the upper- left panel; ^ < 0.1 with thickness of 
15/i^^Mpc in the upper-right panel; z < 0.2 with thickness of 20/i~^Mpc in the 
lower panel. 



^http: / /www.sdss.org/ 
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Figure 22: 3D redshift-space map centered on us, and its projection on the celes- 
tial sphere of SDSS galaxy subset, including the three main regions, (rcf. [31]) 
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Figure 23: Redshift slices of SDSS galaxy data around the equatorial plane. 
The redshift limits and the thickness of the planes are: Upper-left z < 0.05, 
lO/i-^Mpc; Upper-right z < 0.1, 15h~^Mpc; Louver z < 0.2, 20hr^Mpc. The 
size of points has been adjusted. Note that the data for the Southern part are 
sparser than those for the Northern part, especially for thick slices, (rcf. [31]) 
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6.1.3 The 6dF galaxy redshift survey 

The 6dF (6-degree Field) is a survey of redshifts and peculiar velocities of 
galaxies selected primarily in the Near Infrared from new 2MASS (Two Micron 
All Sky Survey) catalogue^. One goal is to measure redshifts of more than 
170,000 galaxies over nearly the whole Southern sky. Another exciting aim of 
the survey is to measure peculiar velocities (using 2MASS photometry and 6dF 
velocity dispersions) of about 15,000 galaxies out to 150/i~"'^Mpc. The high 
quality data of this survey could revive peculiar velocities as a cosmological 
probe (which was very popular about 10-15 years ago). Observations have so 
far obtained nearly 40,000 redshifts and completion is expected in 2005. 

6.1.4 The DEEP galaxy redshift survey 

The DEEP survey is a two-phased project using the Keck telescopes to study the 
properties and distribution of high redshift galaxies ^ . Phase 1 used the LRIS 
spectrograph to study a sample of ~ 1000 galaxies to a limit of 1=24.5. Phase 2 
of the DEEP project will use the new DEIMOS spectrograph to obtain spectra 
of ~ 65,000 faint galaxies with redshifts z ^ 1. The scientific goals arc to 
study the evolution of properties of galaxies and the evolution of the clustering 
of galaxies compared to samples at low redshift. The survey is designed to 
have the fidelity of local redshift surveys such as the LCRS survey, and to be 
complementary to ongoing large redshift surveys such as the SDSS project and 
the 2dF survey The DEIMOS/DEEP or DEEP2 survey will be executed with 
resolution R 4000, and we therefore expect to measure lincwidths and rotation 
curves for a substantial fraction of the target galaxies. DEEP2 will thus also be 
complementary to the VLT/VIRMOS project, which will survey more galaxies 
in a larger region of the sky but with much lower spectral resolution and with 
fewer objects at high redshift. 

6.1.5 The VIRMOS galaxy redshift survey 

The on-going Franco-Italian VIRMOS project^ has delivered the VIMOS spec- 
trograph for the European Southern Observatory Very Large Telescope (ESO- 
VLT). VIMOS is a Visible imaging Multi-Object Spectrograph with outstand- 
ing multiplex capabilities: with 10 arcsec slits, spectra can be taken of 600 
objects simultaneously. In integral field mode, a 6400-flbre Integral Field Unit 
(IFU) provides spectroscopy for all objects covering a 54x54 arcsec^ area. VI- 
MOS therefore provides unsurpassed efficiency for large surveys. The; VIRMOS 
project consists of: Construction of VIMOS, and a Mask Manufacturing Unit 
for the ESO-VLT. The VIRMOS- VLT Deep Survey (VVDS), a comprehensive 

^http://www.mso.anu.edu.au/6dFGS/ 
''http://www.ipac.caltech.edu/2mass/ 
^http://deep.berkeley.edu/ marc/deep/ 
^ http: // www. astrsp-mrs . fr / virmos / 
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imaging and redshift survey of the deep universe based on more than 150,000 
redshifts in four 4 sq.-degree fields. 

6.2 Cosmological parameters from 2dFGRS 

6.2.1 The Power spectrum of 2dF GalELxies on leirge scales 

An initial estimate of the convolved, rcdshift-spacc power spectnim of the 2dF- 
GRS was determined by Percival et al. (2002) [72] for a sample of 160,000 
redshifts. On scales 0.02/i/Mpc < k < 0.15/i/Mpc, the data are fairly robust 
and the shape of the power spectrum is not significantly affected by redshift- 
space distortion or non- linear effects, while its overall amplitude is increased 
due to the linear rcdshift-space distortion effect (see §5). 

If one fits the A-CDM model predictions to the 2dFGRS power spectrum 
(Fig. 24) over the above range in fc, one can constrain the cosmological pa- 
rameters. For instance, assuming a Gaussian prior on the Hubble constant 
h = 0.7 ± 0.07 (from ref. [22]), Percival et al. (2002) [72] obtained the 68 per- 
cent confidence limits on the shape parameter f2„/i = 0.20±0.03, and the baryon 
fraction Ot/Jl,,, = 0.15 ± 0.07. For a fixed set of cosmological parameters, i.e., 
n = 1, f2„ = 1 — = 0.3, f2b^^ = 0.02 and h = 0.70, the rms mass fluctuation 
amplitude of 2dFGRS galaxies smoothed over a top-hat radius of 8/i~^Mpc in 
redshift space turned out to be a§g{Ls, Zs) « 0.94. 

6.2.2 An upper limit on neutrino masses 

The recent results of atmospheric and solar neutrino oscillations [24, 1] imply 
non-zero mass-squared differences of the three neutrino flavours. While these 
oscillation experiments do not directly determine the absolute neutrino masses, 
a simple assumption of the neutrino mass hierarchy suggests a lower limit on the 
neutrino mass density parameter, Cli, = m^^toth~'^ / {QieV) « 0.001. Large scale 
structure data can put an upper limit on the ratio due to the neutrino 

'free streaming' effect [33]. By comparing the 2dF galaxy power-spectrum of 
fluctuations with a four-component model (baryons, cold dark matter, a cos- 
mological constant and massive neutrinos) it was estimated that < 0.13 
(95% CL), or with concordance prior of ri„ = 0.3, fiy < 0.04, or an upper Hmit 
of 2 eV on the total neutrino mass, assuming a prior of /i w 0.7 [17, 18] 
(see Fig. 24). In order to minimize systematic effects due to biasing and non- 
linear growth, the analysis was restricted to the range 0.02 < k < 0.15 /iMpc~^. 
Additional cosmological data sets bring down this upper limit by a factor of 
two [79]. 

6.2.3 Combining 2dFGRS & CMB 

While the CMB probes the fluctuations in matter, the galaxy redshift surveys 
measure the perturbations in the light distribution of particular tracer (e.g. 
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Figure 24: The power spectrum of the 2dFGRS. The points with error bars 
show the measured 2dFGRS power spectrum measurements in redshift space, 
convolved with the window function. Also plotted are linear CDM models with 
neutrino contribution of Cli, = 0, 0.01 and 0.05 (bottom to top lines) The other 
parameters are fixed to the concordance model. The good fit of the linear 
theory power spectrum at fc > 0.15 /iMpc^^ is due to a conspiracy between the 
non-linear gravitational growth and the Finger-of-God smearing, (refs. [72, 17]) 



galaxies of certain type). Therefore, for a fixed set of cosmological parameters, 
a combination of the two can constrain better cosmological parameters, and it 
can also provide important information on the way galaxies are 'biased' relative 
to the mass fluctuations, 

The CMB fluctuations are commonly represented by the spherical harmonics 
Cg. The connection between the harmonic £ and k is roughly 

for a spatially-flat universe. For ri„ = 0.3 the 2dFGRS range 0.02 < k < 
0.15 /iMpc^^ corresponds approximately to 200 < £ < 1500, which is well cov- 
ered by the recent CMB experiments. 

Recent CMB measurements have been used in combination with the 2dF 
power spectrum. Efstathiou et al. (2002) [19] showed that 2dFGRS-FCMB 
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provide evidence for a positive Cosmological Constant fl^ ~ 0.7 (assuming w = 
— 1), independently of the studies of Supernovae la. As explained in ref. [72], 
the shapes of the CMB and the 2dFGRS power spectra are insensitive to Dark 
Energy. The main important effect of the Dark Energy is to alter the angular 
diameter distance to the last scattering, and thus the position of the first acoustic 
peak. Indeed, the latest result from a combination of WMAP with 2dFGRS and 
other probes gives h = 0.7ll° °^, n^,h^ = 0.0224 ± 0.0009, n^h^ = 0.1351:°;°°^, 
as = 0.84 ± 0.04, Qtot = 1-02 ± 0.02, and w < -0.78 (95% CL, assuming 
w > -1) [79]. 

6.2.4 Redshift-space distortion 

An independent measurement of cosmological parameters on the basis of 2dF- 
GRS comes from redshift-space distortions on scales ^ 10h~^Mpc; a correla- 
tion function ^(7r,cr) in parallel and transverse pair separations tt and a. As 
described in §5, the distortion pattern is a combination of the coherent in- 
fall, parameterized by /? = and random motions modelled by an expo- 

nential velocity distribution function (eq.[133]). This methodology has been 
applied by many authors. For instance. Peacock et al. (2001) [66] derived 
P{L, = 0.17, Zs = 1.9L*) = 0.43 ± 0.07 and Hawkins et al. (2003) [30] obtained 
P{Ls = 0.15, Zs = 1.4L*) = 0.49 ± 0.09 and velocity dispersion cTp = 506 ± 52 
km/sec. Using the full 2dF+CMB likelihood function on the (6, n„) plane, La- 
hav et al. (2002) [42] derived a slightly larger (but consistent within the quoted 
error-bars) value, 13(1^ = 0.17, z^ = 1.9L*) ~ 0.48 ± 0.06. 

6.2.5 The bi-spectrum and higher moments 

It is well established that important information on the non-linear growth of 
structure is encoded at the high order moments, e.g. the skewness or its Fourier 
version, the bi-spectrum. Verde et al. [92]) computed the bi-spectrum of 
2dFGRS and used it to measure the bias parameter of the galaxies. They 
assumed a specific quadratic biasing model: 

5g = hSm + 2^2^m- 

By analysing 80 million triangle configurations in the wavenumber range 0.1 < 
k < 0.5VMpc they found bi = 1.04 ± 0.11 and 63 = -0.054 ± 0.08, in sup- 
port of no biasing on large scale. This is a non-trivial result, as the analysis 
covers non-linear scales. Baugh et al. [3] and Croton et al. [14]) measured the 
moments of the galaxy count probability distribution function in 2dFGRG up 
to order p = 6 (order p = 2 is the variance, p = 3 is the skewness, etc.). They 
demonstrated the hierarchical scaling of the averaged p-point galaxy correlation 
functions. However, they found that the higher moments are strongly affected 
by the presence of two massive superclusters in the 2dFGRS volume. This poses 
the question of whether 2dFGRS is a 'fair sample' for high order moments. 
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6.3 Luminosity and spectral-type dependence of galaxy 
clustering 

Although biasing was commonly neglected until the early 1980s, it has become 
evidc;nt observationally that on scales ^ lOhr^Mpc different galaxy popula- 
tions exhibit different clustering amplitudes, the so-called morphology-density 
relation [16]. As discussed in §4, galaxy biasing is naturally predicted from a 
variety of theoretical considerations as well as direct numerical simulations [37, 
59, 15, 86, 87, 91]. Thus in this subsection, we summarize the extent to which 
the galax;y clustering is dependent on the luminosity, spectral-type and color of 
the galaxy sample from the 2dFGRS and SDSS. 

6.3.1 2dFGRS: Clustering per luminosity and spectral type 

Madgwick et al. (2003) [45] applied the Principal Component Analysis to com- 
press each galaxy spectrum into one quantity, r/ w 0.5 pci + pc2. Qualitatively, 
T] is an indicator of the ratio of the present to the past star formation activ- 
ity of each galaxy. This allows one to divide the 2dFGRS into rj-types, and 
to study e.g. luminosity functions and clustering per type. Norberg et al. 
(2002) [61] showed that, at all luminosities, early-type galaxies have a higher 
bias than late-type galaxies, and that the biasing parameter, defined here as 
the ratio of the galaxy to matter correlation function, b = \Jig/£,m varies as 
6/6* = 0.85 -h 0.15I//I/*. Figure 25 indicates that for L* galaxies, the real space 
correlation function amplitude of r} early-type galaxies is ~ 50% higher than 
that of late- type galaxies. 

Figure 26 shows the redshift-space correlation function in terms of the line-of- 
sight and perpendicular to the line-of-sight separation ^(cr, tt). The correlation 
function calculated from the most passively ('red', for which the present rate 
of star formation is less than 10 % of its past averaged value) and actively 
('blue') star-forming galaxies. The clustering properties of the two samples are 
clearly distinct on scales ^ 10 Mpc. The 'red' galaxies display a prominent 
'finger-of-god' effect and also have a higher overall normalization than the 'blue' 
galaxies. This is a manifestation of the well-known morphology-density relation. 
By fitting C(7r, cr) over the separation range 8 — 20/i~^Mpc for each class, it 
was found that fiacUve = 0.49 ± 0.13, (3passive = 0.48 ± 0.14 and corresponding 
pairwise velocity dispersions cTp of 416±76 km/sec and 612±92 km/sec [45]. At 
small separations, the real space clustering of passive galaxies is stronger than 
that of active galaxies, the slopes 7 are respectively 1.93 and 1.50 (Fig. 27) and 
the relative bas between the two classes is a declining function of separation. 
On scales larger than 10 ft," ^ Mpc the biasing ratio is approaching unity. 

Another statistic was applied recently by Wild et al. [94] and Conway et 
al.[13], of a joint counts- in-cells on 2dFGRS galaxies, classified by both color 
and spectral type. Exact linear bias is ruled out on all scales. The counts 
are better fitted to a bivariate lognormal distribution. On small scales there is 
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Figure 25: The variation of the galaxy biasing parameter with himinosity, rel- 
ative to an galaxy for the full sample and for subsamples of early and late 
spectra types (ref. [61]). 



evidence for stochasticity. Further investigation of galaxy formation models is 
required to understand the origin of the stochasticity. 
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Figure 26: The two point correlation function ^(cr, tt) plotted for passively (left) 
and actively (right) star-forming galaxies. The line contours levels show the 
best-fitting model (ref. [45]). 
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Figure 27: The correlation function for early and late spectral types. The solid 
lines show best-fitting models, whereas the dashes lines are extrapolations of 
these lines (ref. [45]). 
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6.3.2 SDSS: Two-point correlation functions per luminosity and color 



Zehavi et al. (2002) [98] analyzed the Early Data Release (EDR) sample of 
the SDSS 30,000 galaxies to explore the clustering of per luminosity and color. 
The inferred real-space correlation function is well described by a single power- 
law; = (r/6.1 ± 0.2/i-iMpc)-i '^5±o.o3^ f^j. Q.l/i-iMpc < r < le/i-^Mpc. 
The galaxy pairwise velocity dispersion is 1T12 ~ 600 ± 100km/ sec for projected 
separations 0.15/i~^Mpc < Vp < 5/i~^Mpc. When divided by color, the red 
galaxies exhibit a stronger and steeper real-space correlation function and a 
higher pairwise velocity dispersion than do the blue galaxies. In agreement with 
2dFGRS there is clear evidence for scale-independent luminosity bias at r ~ 
10/i~^Mpc. Subsamples with absolute magnitude ranges centered on M, — 1.5, 
M^,, and M, -1-1.5 have real-space correlation functions that are parallel power 
laws of slope « —1.8 with correlation lengths of approximately 7.4ft,^^Mpc, 
6.3/i~^Mpc, and 4.7ft.~^Mpc, respectively. 

Figures 27 and 28 pose an interesting challenge to the theory of galaxy 
formation, to explain why the correlation functions per luminosity bins have 
similar slope, while the slope for early type galaxies is steeper than for late 
type. 




(h ' Mpc) 



-23.0 < M < -21.5 
0-100 < z < 174 

-21.5 < M < -20.0 
0.052 < z < 0,097 

-20.0 < M < -18.5 
0.037 < z < 0.051 



(h"' Mpc) 



Figure 28: Upper: The SDSS (EDR) projected correlation function for blue 
(squares), red (triangles) and the full sample, with best-fitting models over the 
range 0.1 < rp < 16/i~^Mpc. Lower: The SDSS (EDR) projected correlation 
function for three volume- limited samples, with absolute magnitude and redshift 
ranges as indicated and best-fitting power-law models, (ref [98]) 
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6.3.3 SDSS: Three-point correlation functions and the nonlinear bi- 
asing of galaxies per luminosity and color 

Let us move next to the the three-point correlation functions (3PCF) of galaxies, 

which arc the lowest-order unambiguous statistic to characterize non-Gaussianities 
due to nonlinear gravitational evolution of dark matter density fields, formation 
of luminous galaxies and their subsequent evolution. The determination of the 
3PCF of galaxies was pioneered by Peebles & Groth (1975) and Groth & Peebles 
(1977) [70, 27] using the Lick and Zwicky angular catalogs of galaxies. They 
found that the 3PCF C(^'12,''23,t'3i) obeys the hierarchical relation: 

C(r-i2,r23,r3i) = 0rK(ri2)e(r-23)+C(^23)^(r3i)+e(r3i)e('^i2)] (170) 

with Qi being a constant. The value of Qi in real space de-projected from 
these angular catalogues is 1.29 ± 0.21 for r < 3/i~^Mpc. Subsequent analyses 
of redshift catalogs confirmed the hierarchical relation, at least approximately, 
but the value of Qz (in redshift space) appears to be smaller, 0.5 ~ 1. 

As we have seen in §6.3.2, galaxy clustering is sensitive to the intrinsic prop- 
erties of the galaxy samples under consideration, including their morphological 
types, colors, and luminosities. Nevertheless the previous analyses were not 
able to examine those dependences of 3PCFs because of the limited number of 
galaxies. Indeed Kayo et al. (2004) [39] are the first to perform the detailed 
analysis of 3PCFs explicitly taking account of the morphology, color and lu- 
minosity dependence. They constructed volume-limited samples from a subset 
of the SDSS galaxy redshift data, 'Large scale Structure Sample 12'. Specif- 
ically they divided each volume limited sample into color subsamples of red 
(blue) galaxies, which consist of 7949 (8329), 8930 (8155), and 3706 (3829) 
galaxies for -22 < M,. - 5 log ft, < -21, -21 < - 51og/i < -20, and 
-20 < Mr - 5 log /i < -19, respectively. 

Figure 29 indicates the dimensionless amplitude of the 3PCFs of SDSS galax- 
ies in redshift space: 

for the equilateral triplets of galaxies. The overall conclusion is that is al- 
most scale-independent and ranges between 0.5 and 1.0, and that no systematic 
dependence is noticeable on luminosity and color. This implies that the 3PCF 
itself does depend on the galaxy properties since two-point correlation functions 
(2PCFs) exhibit clear dependence on luminosity and color. Previous simulations 
and theoretical models [82, 53, 49, 85] indicate that Q decreases with scale in 
both real and redshift spaces. This trend is not seen in the observational results. 
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Figure 29: Dimensionless amplitude of the three-point correlation functions 
of SDSS galaxies in redshift space. The galaxies are classified according to 
their colors; all galaxies in open circles, red galaxies in solid triangles, and blue 
galaxies in crosses, (ref. [39]). 
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Figure 30: Same as Fig. 29 but for the inverse of the biasing parameter defined 
through the two-point correlation functions, (ref. [39]). 
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In order to demonstrate the expected dependence in the current samples, 
they compute the biasing parameters estimated from the 2PCFs: 



b.M^x i ^"'^'\,y (172) 

where the index i runs over each sample of galaxies with different colors and 
luminosities. The predictions of the mass 2PCFs in redshift space, Cz,acdm(s), 
in the A cold dark matter model arc computed following rcf. [28]. 

As an illustrative example, consider a simple bias model in which the galaxy 
density field i5g,i for the i-th population of galaxies is given by 

^S,i = ^g,i(l)'^mass + ^g,i(2) <5mass ■ (1^3) 

If both feg,,;(i) and 6g^j(i) are constant and mass density field Smass ^ 1, equation 
(171) implies that 

Qg,i ~ 7 Qmass ^" 72" • (-^'^^) 

Thus the linear bias model (&g,i(2) = 0) simply implies that (5g,i is inversely 
proportional to ?>g,i(i), which is plotted in Figure 30. Comparison of Figures 
29 and 30 indicates that the biasing in the 3PCFs seems to compensate the 
difference of Qg purely due to that in the 2PCFs. 

Such behavior is unlikely to be explained by any simple model inspired by 
the pcrturbativc expansion like equation (173). Rather it indeed points to a kind 
of regularity or universality of the clustering hierarchy behind galaxy formation 
and evolution processes. Thus the galaxy biasing seems much more complex 
than the simple deterministic and linear model. More precise measurements 
of 3PCFs and even higher-order statistics with future SDSS datasets would be 
indeed valuable to gain more specific insights into the empirical biasing model. 



6.4 Topology of the universe: analysis of SDSS galcixies 
in terms of Minkowski functionals 

All the observational results presented in the preceding subsections were re- 
stricted to the two-point statistics. As emphasized in §3, the clustering pattern 
of galaxies has much richer content than the two-point statistics can probe. 
Historically the primary goal of the topological analysis of galaxy catalogues: 
was to test Gaussianity of the primordial density fluctuations. Although the 
major role for that goal has been superseded by the CMB map analysis [41], 
the proper characterization of the morphology of large-scale structure beyond 
the two-point statistics is of fundamental importance in cosmology. In order to 
illustrate a possibility to explore the topology of the universe by utilizing the 
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new large surveys, we summarize the results of the Minkowski Functionals (MF) 
analysis of SDSS galaxy data [31]. 

In an apparent-magnitude limited catalogue of galaxies, the average num- 
ber density of galaxies decreases with distance because only increasingly bright 
galaxies are included in the sample at larger distance. With the large redshift 
surveys it is possible to avoid this systematic change in both density and galaxy 
luminosity by construct volume-limited samples of galaxies, with cuts on both 
absolute-magnitude and redshift. This is in particular useful for analyses such 
as MF and was carried out in the analysis shown here. 



Figure 31 shows the MFs as a function of Vf defined from the volume frac- 
tion [26]: 



This is intended to map the threshold so that the volume fraction on the high- 
density side of the isodensity surface is identical to the volume in regions with 
density contrast 5 = Vfa, for a Gaussian random field with r.m.s. density fluctu- 
ations a. If the evolved density field may approximately have a good one-to-one 
correspondence with the initial random-Gaussian field, then this transformation 
removes the effect of evolution of the PDF of the density field. Under this as- 
sumption, the MFs as a function of volume fraction would be sensitive only to 
the topology of the isodensity contours rather than evolution with time of the 
density threshold assigned to a contour. While the limitations of the approxi- 
mation of monotonicity in the relation between initial and evolved density fields 
are well recognized [40], we plot the result in this way for simplicity. 

The good match between the observed MFs and the mock predictions based 
on the LCDM model with the initial random-Gaussianity, as illustrated in 
Figure 31, might be interpreted to imply that the primordial Gaussianity is 
confirmed. A more conservative^ interpretation is that, given the size of the 
estimated uncertainties, these data do not provide evidence for initial non- 
Gaussianity, i.e., the data are consistent with primordial Gaussianity. Unfortu- 
nately due to the statistical limitation of the current SDSS data, it is not easy to 
put more quantitative statement concerning the initial Gaussianity. Moreover, 
in order to go further and place more quantitative constraints on primordial 
Gaussianity with upcoming data, one needs a more precise and reliable theo- 
retical model for the MFs which properly describes the nonlinear gravitational 
effect possibly as well as galaxy biasing beyond the simple mapping on the 
basis of the volume fraction. In fact, galaxy biasing is a major source of uncer- 
tainty for relating the observed MFs to those obtained from the mock samples 
for dark matter distributions. If LCDM is the correct cosmological model, the 
good match of the MFs for mock samples from the LCDM simulations to the 
observed SDSS MFs may indicate that nonlinearity in the galaxy biasing is rel- 
atively small, at least small enough that it does not significantly affect the MFs 
(the MFs as a function of fa remain unchanged for the linear biasing). 




(175) 
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Figure 31: MFs as a function of vi for Rq = 5h~^Mpc for SDSS data. Averaged 
MPs of the mock samples are plotted for LCDM (solid lines) and SCDM (long 
dashed lines). Gaussian model predictions (eqs.[112] to [115]) are also plotted 
with short dashed lines. The results favor the LCDM model with random- 
Gaussian initial conditions, (ref. [31]). 
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6.5 Other statistical measures 

In this section, we have presented the results on the basis of particular statis- 
tical measures including the two-point correlation functions, power spectrum, 
rcdshift distortion, and Minkowski functionals. Of course there arc other useful 
approaches in analysing redshift surveys; the void probability function, counts- 
in-cells, Voronoi cells, percolation, and minimal spanning trees. Another area 
not covered here is optimal reconstruction of density field (e.g. using the Wiener 
filtering). The reader is referred to a good summary of those and other method- 
ology in the book by Martinez & Saar (2002) [48]. 

Admittedly the results that wc presented here arc rather observational and 
phenomenological, and far from being well-understood theoretically. It is quite 
likely that when other on-going and future surveys are be analysed in great 
detail, the nature of galaxy clustering will be revealed in a much more quan- 
titative manner. They are supposed to act as a bridge between cosmological 
framework and galaxy formation operating in the universe. While the proper 
understanding of physics of galaxy formation is still far away, the future redshift 
survey data will present interesting challenges for constructing models of galaxy 
formation. 
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7 Discussion 



As a classical probe, galaxy redshift surveys still remain an important tool for 
studying cosmology and galaxy formation. On the large scales (> Wh~^ Mpc or 
so) they nicely complement the Cosmic Microwave Background, Supernovae la 
and Gravitational Lensing in quantifying in detail the cosmological model. On 
the small scales (< 10h~^ Mpc) the clustering patterns of different galaxy types 
(defined by structural or spectral properties) provide important constraints on 
models of biased galaxy formation. 

The redshift surveys mainly constrain fl^ via both redshift distortion (which 
also depends on biasing) and the shape of the A-CDM power spectrum, which 
depends on the primordial spectrum, the product Om/i and also the baryon 
density fib- Redshift surveys at a given epoch are not sensitive to the Dark 
Energy (or the Cosmological Constant, in a specific case), but combined with 
the CMB they can constrain the cosmic Equation of State. 

A good example of the importance of the redshift surveys in Cosmology is 
given by the recent WMAP analysis of cosmological parameters, where the es- 
timation of certain parameters was much improved by adding the 2dF power 
spectrum of fluctuations ([79]) or the SDSS power spectrum ([88]). This is il- 
lustrated in Table 2 by contrasting the WMAP alone derived parameters from 
those derived from WMAP + SDSS ([88]). Such results are sensitive to the as- 
sumed parameter space and priors, but for simplicity we quote here the results 
for the simple six-parameter model. In this analysis it was assumed that the 
universe is flat, the fluctuations are adiabatic, no gravity waves, no running tilt 
of the spectral index, negligible neutrino masses and that the Dark Energy is 
in the form of Einstein's Cosmological constant {w = —1). Within the A-CDM 
model this scenario can be characterised by the six parameters given in Table 2 
based on [88]. The WMAP data used are both the temperature and polarization 
fluctuations. It can bee seen that by adding the SDSS information more than 
halves the WMAP-only error bars on some of the parameters. These results are 
in good agreement with the joint analysis WMAP-|-2dF ([79]). 

We emphasize that these parameters were fitted assuming the A-CDM model. 
While the degree of such phenomenological successes of the A-CDM model is 
truly amazing, there are many fundamental open questions: 

(i) Both components of the model, A and CDM, have not been directly mea- 
sured. Are they 'real' entities or just 'epicycles'? Historically epicycles 

were actually quite useful in forcing observers to improve their measure- 
ments and theoreticians to think about better models! 

(ii) 'The Old Cosmological Constant problem': Why is at present so small 
relative to what is expected from Early Universe physics? 

(iii) 'The New Cosmological Constant problem': Why is fl^ ~ il\ at the 
present-epoch ? Why is w — 1 ? Do we need to introduce new physics 
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Table 2: Derived 6 cosmological parameters from WMAP alone (temperature 
& polarization) vs. WMAP+SDSS, from Tables 2 and 3 of Tegmark et al. 
(2003). [88] The quoted error bars correspond to l-tr. 



Symbol 


WMAP alone 


WMAP+SDSS 


description 


Qa 

<78 
Us 
T 


'O-O.IO 
U.U/40_o 0019 
U.i4U_o 018 

0.991°;}! 

i-02l°;J^ 

0.2ll°;?t 


699+" °*^ 

U.U^O^_o,ooio 
U.14O4_0 0082 

n qi 7+0.090 
n q77+0.039 

' -0.025 
r. 124+0.083 


dimensionless cosmological constant 
baryon density parameter 
total matter density parameter 
mass fluctuation amplitude 
at 8/i~^Mpc sphere (linear theory) 
primordial scalar spectral index 
at = 0.05/Mpc 
re-ionization optical depth 



or to invoke the Anthropic Principle to explain it ? 

(iv) There are still open problems in A-CDM on the small scales e.g. galaxy 

profiles and satellites. 

(v) Could other (yet unknown) models fit the data equally well ? 

(vi) Where does the field go from here ? Should the activity focus on refinement 
of the cosmological parameters within A-CDM, or on introducing entirely 
new paradigms ? 

Even if ACDM turns out to be the correct model, it is not yet the "end 
of cosmology". Beyond the 'zero-th order' task of finding the cosmological 
parameters of the FRW model, we would like to understand the non-linear 
growth of mass density fluctuations and then the formation and evolution of 
luminous objects. The wealth of data of galaxy images and spectra in the new 
surveys calls for the development of more detailed models of galaxy formation. 
This is important so the comparison of the measurements (e.g. correlation 
function per spectral type or colour) and the models could be done on equal 
footing, with the goal of constraining scenarios of galaxy formation. There is 
also room for new statistical methods to quantify the 'cosmic web' of filaments, 
clusters of voids, for effective comparison with the simulations. It may well be 
that in the future the cosmological parameters will be fixed by the CMB, SN la 
and other probes. Then, for fixed cosmological parameters, one may use redshift 
surveys primarily to study galaxy biasing and evolution with cosmic epoch. 
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